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In  this  report,  the  reflection  and  scattering  models  and  the  propagation  model 
including  the  effects  of  realistic  ocean  boundaries  are  further  developed.  ^U-I 


A  brief  introduction  which  includes  a  short  history  of  the  scattering  and 
propagation  work  performed  at  ARL,  and  the  considerations  which  motivated  the 
research  represented  in  this  report  is  given  in  Chapter  I.  In  Chapter  II, 
the  theoretical  foundations  of  the  scattering  by  acoustically  penetrable 
surfaces  are  presented.  The  formulas  developed  in  Chapter  II  are  used  in 
Chapter  III  to  investigate  the  effects  of  various  realistic  ocean  bottoms  on 
the  performance  of  bottom  bounce  sonars.  Once  the  effect  of  a  known  bottom 
on  the  signal  can  be  predicted,  the  problem  of  p:  edicting  the  physical  pro¬ 
perties  of  the  bottom  from  the  signals  it  returns  becomes  manageable. 

Certain  relevant  aspects  of  this  problem  are  discussed  Chapter  IV,  specifi¬ 
cally,  the  prediction  of  the  forward  reflected  signaJ  from  measurements  of 
the  backscattered  signal.  In  Chapter  V,  propagation  in  a  surface  duct  with 
realistic  boundaries  is  treated,  and  the  attenuation  per  mode  due  to  boundary 
roughness  is  calculated  for  the  Epstein  model.  Finally,  the  preliminary 
work  on  the  experimental  measurements  of  the  scattering  by  acoustically 
penetrable,  rough  surfaces  is  presented  in  Appendix  A.  (U-FOUO) 
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ABSTRACT 


In  this  report,  the  reflection  and  scattering  models 
and  the  propagation  model  including  the  effects  of 
realistic  ocean  boundaries  are  further  developed. 

A  brief  introduction  -which  includes  a  short  history  of 
the  scattering  and  propagation  work  performed  at  AEL 
and  the  considerations  which  motivated  the  research 
represented  ir  this  report  is  given  in  Chapter  I.  In 
Chapter  II,  the  theoretical  foundations  of  the  scattering 
by  acoustically  penetrable  surfaces  are  presented.  The 
formulas  developed  in  Chapter  II  are  used  in  Chapter  III 
to  investigate  the  effects  of  various  realistic  ocean 
bottoms  on  the  performance  of  bottom  bounce  sonars. 

Once  the  effect  of  a  known  bottom  on  the  signal  can  be 
predicted,  the  problem  of  predicting  the  physical  prop¬ 
erties  of  the  bottom  from  the  signals  it  returns 
becomes  manageable.  Certain  relevant  aspects  of  this 
problem  are  discussed  in  Chapter  IV,  specifically,  the 
prediction  of  the  forward  reflected  signal  from  measure¬ 
ments  of  the  backscattered  signal.  In  Chapter  V,  propa¬ 
gation  in  a  surface  duct  with  realistic  boundaries  is 
treated,  and  the  attenuation  per  mode  due  to  boundary 
roughness  is  calculated  for  the  Epstein  model.  Finally, 
the  preliminary  work  on  the  experimental  measurements  of 
the  scattering  by  acoustically  penetrable,  rough  surfaces 
is  presented  in  Appendix  A. 
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I .  INTRODUCTION 


The  detection,  localization,  and  classification  of  targets  in 
the  ocean  environment  require  extensive  and  detailed  knowledge  of  the 
propagation  of  acoustic  energy  in  the  ocean.  Typically,  sonar  opera¬ 
tion  is  a  complex  problem  involving  the  propagation  of  sound  in  an 
inhomogeneous  medium  with  imperfect  boundaries  (a  perfect  boundary  is 
a  plane,  acoustically  impenetrable  surface).  As  a  consequence  of  the 
medium  inhomogeneity  and  the  boundary  imperfections,  two  related 
limitations  on  sonar  operation  arise.  Due  to  the  inhomogeneity  of  the 
ocean,  so  called  "shadow  zones"  and  sound  channels  exist  in  the  ocean. 
The  propagation  of  signals  into  a  shadow  zone  or  in  sound  channels 
depends  in  a  very  complicated  manner  on  the  ocean  surface  and  bottom. 
The  effects  of  the  imperfect  ocean  boundaries  are  seen  as  losses  in 
the  strength  of  the  returned  signal  and  as  reverberation  returns.  The 
bottom  bounce  mode  and  the  surface  duct  mode  are  two  examples  of  sonar 
operational  modes  which  are  required  to  obtain  coverage  of  certain 
regions  of  the  ocean  volume,  but  which,  nevertheless,  are  limited  by 
the  ocean  boundaries. 

The  development  of  more  effective  sonars  requires  increased 
knowledge  of  the  physical  nature  of  the  ocean  surface,  the  bottom,  the 
water  between  these  surfaces,  and  the  effect  of  these  factors  on  the 
propagation  of  sound.  Accordingly,  extensive  full  scale  sea  tests 
have  been  conducted  to  study  the  ocean  environment,  the  acoustic  scat¬ 
tering  properties  of  the  ocean  surface  and  bottom,  and  long  range 
propagation  in  the  ocean. 
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Full  scale  experiments  at  sea  are  an  invaluable  source  of  data 
for  typical  conditions  that  might  be  expected  in  an  operational  situa¬ 
tion.  Environmental  conditions,  surface  wave  structure,  and  bottom 
composition  and  roughness  can  only  be  obtained  by  in  situ  measurement. 
However,  the  generally  uncontrolled  sea  conditions  and  lack  of  precise 
knowledge  of  the  test  geometry  can  introduce  ambiguities  when  particular 
aspects  of  scattering  or  propagation  theory  are  to  be  investigated. 

In  response  to  these  difficulties  measurements  using  models  where  the 
geometry  and  the  physical  composition  and  topography  of  the  boundary 
surfaces  could  be  precisely  controlled  were  initiated. 

In  view  of  the  complex  and  entangled  effects  of  the  medium  and 
its  boundaries,  the  treatment  of  acoustic  propagation  in  the  ocean 
is  greatly  facilitated  by  separating  it  into  two  problems.  The  first 
problem  is  the  study  of  propagation  in  an  inhomogeneous  medium  with 
perfectly  reflecting  plane  boundaries.  The  second  problem  deals  with 
propagation  in  a  homogeneous  medium  with  imperfect  boundaries,  that  is, 
the  reflection  and  scattering  problem.  The  reflection  and  scattering 
problem  is  further  simplified  by  separating  the  effects  due  to  the 
topography  of  the  boundary  and  the  effects  due  to  the  acoustic  pene¬ 
trability  of  the  boundary,  which  depends  on  the  physical  composition  of 
the  boundary  media. 

In  1961,  a  research  program  was  initiated  at  Applied  Research 
Laboratories  to  investigate  experimentally  and  theoretically  the  reflec¬ 
tion  and  scattering  of  sound  from  the  ocean  boundaries.  The  first 
studies  conducted  under  this  program  involved  the  reflection  of  sound 
by  plane  layered  sediments  of  the  type  usually  found  in  abyssal  plain 
areas  and  the  scattering  of  sound  by  a  pressure  release  sinusoidal 
surface.  These  studies  were  performed  under  Contract  NObsr-72627  and 
resulted  in  very  satisfactory  agreement  between  the  experimental  and 
theoretical  results. 
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Upon  completion  of  these  investigations,  the  study  of  the 
scattering  properties  of  pressure  release  randomly  rough  surfaces  -was 
undertaken.  This  work  was  done  under  Contracts  NObsr-95127, 
N00024-68-C-1112,  and  N00024-68-C-1275  and  covered  the  period  from 
January  1967  to  April  1970*  The  experimental  program  resulted  in  the 
measurement  of  the  forward,  specular,  and  backscattering  properties  at 
the  frequencies  100,  200,  and  500  kHz  of  four  rough  surfaces  with 
different  rms  heights.  The  frequency-roughness  regimes  investigated 
represent  the  extremes  of  the  cases  where  the  wavelength  is  much  greater 
than  the  roughness  to  the  case  wl  re  the  wavelength  is  very  much  less 
than  the  roughness. 

The  theoretical  investigations  conducted  under  this  program  have 
resulted  in  a  totally  new  and  very  powerful  formulation  of  scattering 
problems.  The  theory  which  is  now  available  gives  very  good  agreement 
for  all  cases  of  forward  and  specular  scattering  and  predicts  the  way 
in  which  the  range  behavior  of  scattering  depends  on  roughness.  The 
theory  had  previously  proved  inadequate  for  backscattering  predictions 
at  low  grazing  angles.  This  defect  in  the  theory  has  now  been  corrected. 

This  final  report  is  a  summary  of  the  research  on  acoustically 
penetrable  rough  surface  scattering  and  propagation  in  a  surface  duct 
with  a  realistically  rough  boundary.  This  work  was  performed  under 
Contract  N00024-70-C-1279  from  1  April  1970  to  51  March  1971.  The 
work  done  under  this  contract  represents  an  effort  to  combine  the  topo¬ 
graphical  and  physiological  effects  of  the  bottom  into  a  unified  scat¬ 
tering  treatment.  Further,  the  surface  duct  propagation  study  is  an 
attempt  to  incorporate  a  realistically  rough  boundary  into  the  propa¬ 
gation  model. 

This  final  report  is  composed  of  the  following  sections.  A 
brief  introduction  to  the  report  has  been  given  in  Chapter  I.  The 
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theoretical  foundations  of  scattering  by  penetrable  surfaces  are 
presented  in  Chapter  II.  In  Chapter  III,  the  formulas  developed  in 
Chapter  II  are  used  to  investigate  the  effects  of  various  realistic 
ocean  bottoms  on  the  operation  of  bottom  bounce  sonars.  Once  the  effect 
of  a  known  bottom  on  the  signal  can  be  predicted,  the  problem  of  pre¬ 
dicting  the  physical  properties  of  the  bottom  from  the  signals  it 
returns  becomes  manageable.  Certain  relevant  aspects  of  this  problem 
are  discussed  in  Chapter  IV,  specifically,  the  prediction  of  the  for¬ 
ward  reflected  signal  from  measurements  of  the  backscattered  signal. 

In  Chapter  V,  the  propagation  in  a  surface  duct  with  realistic  bounda¬ 
ries  is  treated,  and  the  attenuation  per  mode  due  to  boundary  roughness 
is  calculated  for  the  Epstein  model.  Finally,  the  preliminary  work 
on  the  experimental  measurements  of  the  scattering  by  acoustically 
penetrable,  rough  surfaces  is  presented  in  Appendix  A. 

The  theoretical  and  experimental  methods  developed  in  these 
investigations  should  prove  useful  in  the  choice  of  some  of  the  opera¬ 
tional  characteristics  of  present  and  future  sonar  systems.  In 
addition,  further  understanding  of  the  limitations  imposed  by  the 
properties  of  the  boundaries  on  detection  and  localization  capabilities 
of  modern  sonars  should  result  from  these  studies.  Another  important 
consequence  of  these  investigations  is  the  verification  or  improvement 
of  mathematical  models  being  used  in  numerous  computer  simulation 
studies  of  optimum  sonar  operation. 
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II.  REFLECTION  COEFFICIENTS  FOR 
A  ROUGH,  ACOUSTICALLY  PENETRABLE  SURFACE 

A.  Introduction 

It  was  pointed  out  in  Chapter  I  that  the  scattering  Behavior  of 
the  ocean  boundaries  depends  on  both  the  topography  and  the  physical 
composition  of  the  boundary.  The  topography  controls  the  scattering 
while  the  physical  composition  of  the  boundaries  determines  the 
acoustic  penetrability  of  a  surface,  and  as  will  be  seen  later,  the 
two  effects  interact.  On  some  boundaries,  for  example  the  water-air 
interface,  the  acoustic  impedance  is  so  large  that  practically  all  of 
the  sound  remains  trapped  in  the  ocean;  consequently,  the  scattering 
depends  entirely  on  the  surface  topography  (including  bubbles,  etc.). 

On  the  other  hand,  the  problem  of  scattering  at  the  ocean  bottom  requires 
for  its  solution  the  inclusion  of  the  effects  of  both  the  topography 
and  the  physical  composition  of  the  bottom.  The  scattering  due  to 
topography  alone  has  a  very  extensive  literature  and  is  now  thought  to 
be  well  understood.  However,  very  little  is  known  about  the  effect  of 
the  acoustic  penetrability  of  the  boundary  on  scattering.  Most  of  the 
literature  available  treats  only  the  single  interface  case  and  does  not 
give  a  very  good  account  of  the  interrelation  of  the  topography  and 
acoustic  penetrability  of  the  scattering  surface.  In  this  section  of 
the  final  report,  the  effect  of  the  acoustic  penetrability  of  the  ocean 
bottom  will  be  incorporated  into  the  scattering  integral  developed 
earlier  in  the  program,  and  the  interaction  of  the  topography  with  the 
penetrability  will  be  discussed. 
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Historically  the  problem  of  calculating  the  scattered  pressure 
from  an  acoustically  penetrable  interface,  such  as  the  ocean  bottom, 
has  been  treated  in  a  number  of  ways.  The  most  common  method  has  been 
to  assume  that  the  pressure  scattered  by  the  acoustically  penetrable 
interface  is  given  by  the  pressure  scattered  from  a  perfectly  reflecting 
surface  modified  by  a  suitable  reflection  coefficient  to  account  for  the 
penetration  losses.  The  form  of  the  reflection  coefficient  is  deter¬ 
mined  by  the  boundary  conditions  at  the  surface,  and  its  value  is  deter¬ 
mined  by  the  local  geometry  (local  slope  of  the  surface).  If  the  surface 
is  plane  or  suitably  smooth,  the  local  geometry  is  the  same  everywhere 
on  the  surface  and  the  reflection  coefficient  is  usually  assumed  to  be 
the  Rayleigh  coefficient.  It  should  be  pointed  out  that  to  arrive  at 
the  Rayleigh  coefficient  it  is  necessary  to  assume  that  the  incoming 
radiation  is  a  plane  wave  and  that  the  outgoing  (reflected)  wave  is 
radiated  in  the  local  specular  direction.  This  assumption  is  incon¬ 
sistent  with  the  integral  equation  of  scattering  which  is  based  only  on 
wave  mechanics  (as  opposed  to  ray  theory).  It  will  be  shown  that  the 
restriction  on  the  outgoing  wave  can  be  easily  removed  in  the  context 
of  the  potential  formulation,  and  in  the  limit  kr  »  1  (k  is  the  wave 
number  and  r  is  the  separation  distance  of  the  source  and  surface), 
this  formulation  results  in  a  reflection  coefficient  identical  to  the 
Rayleigh  coefficient  for  the  single  rough  interface  case. 

For  smooth  surfaces,  it  is  usually  assumed  that  the  reflection 
coefficient  can  be  removed  from  the  scattered  pressure  integral.  Since 
the  slopes  are  small,  the  local  value  of  the  reflection  coefficient  may 
be  everywhere  approximated  by  the  value  of  Rayleigh  coefficient  refer¬ 
enced  to  the  mean  plane,  and  this  coefficient  can  be  removed  from  the 
integral  as  a  constant.  When  this  is  done,  the  remaining  integral  can 
be  evaluated  in  a  straightforward  manner,  and  it  gives  the  expression 
for  the  field  scattered  by  a  pressure  release  surface.  However,  if 
the  surface  is  rough,  the  boundary  conditions  are  only  given  locally, 
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so  the  value  of  the  reflection  coefficient  varies  from  point  to  point  on 
the  surface.  Thus,  the  reflection  coefficient  cannot  be  removed  from 
the  integral  for  the  scattered  pressure  except  in  an  averaged  form.  Due 
to  the  large  rms  roughness  of  the  surfaces  considered  in  this  program, 
a  treatment  of  the  reflection  coefficient  within  the  context  of  the 
scattered  pressure  integral  has  been  pursued.  This  basic  approach  has 
been  modified  in  several  ways  to  give  expressions  for  the  reflection 
coefficient.  These  results  will  be  discussed  later  in  this  section 
and  may  be  generally  classified  as  the  series,  the  exact,  and  the  sta¬ 
tionary  phase  expressions.  In  addition,  the  reflection  coefficient 
given  by  Kuo'*'  is  included  in  this  section  for  the  sake  of  comparison. 

B.  Formal  Development 

The  reflection  coefficients  presented  in  this  section  were  all 

calculated  for  the  same  sets  of  bottom  parameters.  The  parameters 

2 

used  were  those  given  by  Mackenzie  for  four  separate  coastal  bottom 
samples  (these  data  were  originally  published  by  Hamilton  et  al.^). 

The  samples  used  in  this  section  were:  sample  number  13  (sand-silt- 
clay),  sample  number  27  (clayey  fine  silt),  sample  number  28  (silty 
very  fine  sand),  and  sample  number  32  (sand).  The  attenuations  for 
these  samples  at  30  kHz  were  5*61,  1.42,  2.63,  and  4.48  dB  per  foot, 
respectively. 

1.  Kuo's  Reflection  Coefficient 


The  reflection  coefficient  given  by  E.  7.  T.  Kuo  is  presented 
in  this  section  because  it  is  a  fairly  well  known  result  which  can 
serve  as  a  basis  for  comparison  with  other  reflection  coefficients.1 
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The  coefficients  are  calculated  for  a  backscattering  geometry  on  the 
assumption  of  a  single  interface  and  are  given  by 


»(7) 


P2C27  "  plCin 

P2c27  +  Pi0!11 


+  2 


(i-r'i)p2c|(p2-p1) 

(p2  c27  +  P-jC^) 


(1) 


where 


and 


(2) 


y  =  sin<t> 

Here,  $  is  the  receiver  grazing  angle,  p^  and  c^  are  the  density  and 
sound  velocity  in  the  first  medium  (water),  and  p2  and  are  the 
density  and  sound  velocity  in  the  second  medium  (bottom).  This  coef¬ 
ficient  is  restricted  to  fairly  smooth  surfaces.  Vne  coefficient 
given  by  Eq.  (l)  was  evaluated  using  the  parameters  of  all  four  bottom 
samples,  with  the  results  shown  in  Dwg.  AS-7I-213.  Near  normal  inci¬ 
dence  these  curves  approach  the  value  of  the  Rayleigh  coefficient; 
however,  at  low  grazing  angles  they  become  quite  large.  Although  these 
coefficients,  when  used  in  conjunction  with  Kuo’s  scattering  coeffi¬ 
cients,  give  fairly  good  fits  to  some  backscatter  data,  it  is  thought 
that  reflection  coefficients  that  have  average  magnitudes  greater  than 
unity  over  large  ranges  of  grazing  angle  must  be  regarded  as  nonphysi¬ 
cal.  At  any  given  grazing  angle  a  reflection  coefficient  could,  without 
violating  any  energy  considerations,  have  a  magnitude  greater  than 
i:nity  due  to  focusing  effects.  However,  to  exhibit  this  behavior  over 
a  large  range  of  grazing  angles  implies  a  large  radius  of  curvature  of 
the  mean  plane,  which  is  inconsistent  with  the  assumptions  that  have 
already  been  made  about  the  nature  of  the  surface.  In  conclusion,  to 
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understand  the  relatively  good  agreement  of  Kuo's  results  with 
experimental  data,  it  is  necessary  to  consider  the  entire  expression 
for  the  scattered  field.  The  separation  of  the  scattering  coefficient 
from  the  reflection  coefficient  does  not  appear  to  be  consistent 
because,  ir  the  limit  of  rigid  or  pressure  release  surfaces,  the  reflec¬ 
tion  coefficient  obviously  does  not  give  the  correct  values.  Thus, 

Kuo's  result  does  not  conform  to  the  generally  accepted  definition 
of  a  reflection  coefficient  and  can  only  be  meaningful  in  the  context 
of  his  scattering  theory. 

2.  The  Calculation  of  the  Pressure  Density  in  the  Case  of 
Penetrable  Interfaces 


a.  Single  Interface 

In  the  introduction  to  this  section  it  was  pointed  out 

that  the  adjustment  of  the  boundary  conditions  to  account  for  the 

acoustic  penetrability  of  the  scattering  interface  has  usually  been 

done  by  incorporating  the  Rayleigh  reflection  coefficient  R(cp.  )  into 

1 

the  scattered  pressure 


p  =  R(cp.  )p. 

*s  v  1  *1 

Further,  it  was  stated  that  this  method  is  inconsistent  with  the  wave 
mechanics  viewpoint  of  the  integral  equation.  It  is  a  simple  matter 
to  avoid  this  inconsistency  by  starting  with  the  potential  expression 
for  the  reflected  or  scattered  field  and  then  determining  the  unknown 
pressure  density  on  the  surface. 

Consider  the  geometry  in  Dwg.  AS-70-795*  The  line 
indicated  by  p^  connects  the  source  Q  and  the  surface  point  N.  The 

line  labeled  p  connects  the  surface  point  N  to  the  receiver  at  A. 

R 
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Then  the  incident  pressure  density  at  N  due  to  the  source  at  Q  is 
given  by 


Pi(N)  = 


(3) 


where  r  =  |QN|.  The  reflected  or  scattered  field,  Pr/A),  at  the 
O  K 

receiver  A,  due  to  the  distribution  of  pressure  on  the  surface,  is 
given  by 


pE(A>  s//4  ’  (4) 

s 

where  P  is  some  as  yet  unspecified  pressure  density  (potential)  on 
the  side  of  the  surface  towards  A,  here  r^  =  | AN | ,  k^  is  the  wave 
number  in  the  upper  medium,  and  n  is  the  outward  normal  with  respect 
to  the  surface. 


The  problem  is  to  calculate  the  reflected  or  scattered 

pressure  p_(A)  at  A.  The  calculation  of  p  (a)  could  be  performed  if 
K  K 

P  were  known;  hence,  the  problem  reduces  to  that  of  determining  P. 

Consider  Eq.  (4);  if  the  field  point  A  is  allowed  to 
approach  the  boundary  point  N,  then  Eq.  (4)  becomes 


f  f  /  s\ 

pR(N)  =  2kp(N)  +  J J  P  - -jdS  ,  (5) 

where  the  integral  term  represents  the  multiple  scattering  term. 

If  it  is  assumed  that  as  a  first  approximation  p^N)  s  Rp.(N), 

K  1 
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where  R  is  some  unknown  constant,  then  Eq.  (5)  may  be  solved  by 
iteration  yielding 


prC)  =  V> +  fe  //■*  i<r>tr 


-ik^r  . 
1 


dS 


(6) 


If  multiple  scattering  is  negligible,  then  Eq.  (6)  may  be  approximated 
as 


Pr(N)  s  Rp^(N)  ,  (7) 

and  then 

H  P±(H)  •  (8) 

It  may  be  said  that  Eqs.  (7)  and  (8)  are  fairly  obvious, 
but  in  obtaining  them  in  this  manner,  it  is  clear  that  they  are 
approximations  and  that  the  method  of  performing  the  exact  calculation 
is  indicated. 


The  problem  has  now  been  reduced  to  finding  R,  which 
is  some  pressure  density  coefficient.  It  will  now  be  expedient  to 
collect  certain  quantities  which  will  prove  useful  in  the  calculation 
of  R. 


The  pressure  transmitted  through  the  interface  is  also 
an  outgoing  wave,  and  it  is  given  by 


P 


T 


y 


(9) 


13 


where 


is  taken  at  some  field  point  B  in  the  lower  medium, 
co  is  the  pressure  density  on  the  lower  side  of  the  surface, 
r2  =  |NB|, 

kg  is  the  wave  number  in  medium  2,  and 
n  is  the  outward  normal  with  respect  to  the  surface. 

Using  the  same  analysis  as  in  Eq.  (4),  it  is  easily 

seen  that 


PT  «  Tp^N) 


® s  h  pi(N)  ' 


(10) 


where  T  is  again  some  unknown  constant. 

Now,  for  penetrable,  liquid  boundary  media,  the 
boundary  condition  is  the  continuity  of  the  pressure  across  the  bound¬ 
ary  interface.  A  consequence  of  this  condition  is  Snell's  Law 


k^  *  t  =  kg  •  £ 


(11) 


where  k^  and  kg  are  the  wave  vectors  in  media  1  and  2,  and  £  is  the 
unit  tangent  at  point  N  of  the  boundary.  If  the  interface  were  per¬ 
fectly  plane,  then  Eq.  (ll)  reduces  to 


k  cos©t  =  k„cos9D 

L  L  CL  d 


which  is  the  usual  expression  of  Snell's  Law. 


The  boundary  condition  itself  is  not  sufficient  to 
uniquely  determine  p;  however,  a  further  condition  is  obtained  at  the 
boundary  from  the  hydrodynamic  equation  of  continuity.  This  equation 
implies  that  if  the  interface  remains  intact,  then  the  normal  component 
of  the  velocity  is  continuous  across  the  interface. 


It  is  well  known  that  for  a  homogeneous  medium  the 

velocity  fields  are  easily  related  to  the  pressure  fields.  In  fact, 

the  normal  component  of  the  velocity  field  9 .  in  the  jth  medium  is 

0 

related  to  the  pressure  p^  by  the  simple  expression 


(12) 


where  p .  is  the  static  density  of  jth  medium. 

u 


The  pressure  densities  p  and  o>  are  related  similarly  to 
the  normal  component  of  the  surface  velocity  density  x 


1 

*3  = 

J 


The  application  of  the  boundary  condition  and  the  continuity 
condition  gives  the  two  following  independent  equations: 

p.(N)  +  pR(H)  =  PT(U) 


and 


(13) 


9.(N)  +  9r(N)  =  9t(N) 

Substituting  Eq.  (12)  into  Eqs.  (3),  (^),  and  (9),  respectively,  and 
distinguishing  between  outgoing  and  incoming  waves  give  the  following 
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expressions  for  the  normal  component  of  the  incident,  reflected  and 
transmitted  velocity  fields: 


and 


ik  Sr 

„i(»)  -  — 


5ro  r,  i  ] 

V  L  '  oJ 


Pi 


The  quantity  k^Sr^/Sn  can  be  identified  as  k^  •  n,  the  component  of  the 
wave  vector  in  medium  1  which  is  locally  normal  to  the  surface.  Hence, 


k„i  -\sr-  kisin(ei+fS) 


(15) 


where  P  is  the  inclination  of  the  local  slope  to  the  mean  plane. 
Similarly,  the  term  kgSr^/Sn  can  be  identified  as  the  normal  component 
of  the  wave  vector  in  medium  2.  Further  using  Snell's  Law,  this  com¬ 
ponent  is  uniquely  determined  as 


n2 


=  k. 


Sr2 
2  ^T 


=  kosin(0o+p) 


(16) 


16 


Substituting  Eqs.  (7),  (10),  and  (l4)  into  Eq.  (lj)  and 
solving  for  R  give  the  pressure  density,  n, 


where 


Rp.(N) 
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In  the  limit  as  k,  r  »  1  and  k,r  »  1,  Eq.  (17)  has  the  same  form 
as  the  Rayleigh  reflection  coefficient  for  a  single  interface.  Equa¬ 
tion  (17)  has  certain  inherent  advantages  over  the  Rayleigh  reflection 
coefficient.  First  of  all,  since  the  reflected  or  transmitted  fields 
were  given  in  terms  of  an  integral  equation  which  required  only  the 
incident  pressure  density  on  the  surface,  no  condition  on  the  reflected 
or  transmitted  field  arises.  For  example,  the  condition  of  specularity 
of  the  locally  reflected  field  is  not  present;  this  is  in  accordance 
with  the  usual  wave  mechanics  viewpoint.  Second,  it  was  not  necessary 
to  make  the  appeal  to  ?he  tangent  plane  approximation.  Further,  the 
theory  utilized  spherical  waves  and  contained  no  requirement  that  the 
surface  be  in  the  farfield  of  the  source.  Realistically,  for  most 
sources,  the  surface  should  be  far  enough  removed  so  that  nearfield 
source  behavior  can  be  neglected  in  the  pressure  density  on  the  surface. 


From  the  preceding  derivation,  it  is  clear  that,  for  a 
source  removed  many  wavelengths  from  a  penetrable  surface,  a  suitable 
approximation  for  the  surface  pressxire  potential,  |J,  is  given  by 
modifying  the  incident  pressure  density  p^(N)  by  the  factor  R.  Further 
this  factor  R  is  equivalent  to  the  Rayleigh  reflection  coefficient  but 
does  not  contain  the  explicit  assumption  of  specularity  to  determine 
the  surface  pressure  density  of  the  reflected  field. 


b. 


Double  Interface 


It  was  just  shown  that,  for  a  single  rough  interface, 

the  Rayleigh  reflection  coefficient  referenced  with  respect  to  the  local 

incident  angle  will  be  the  correct  result  if  multiple  scattering  among 

features  on  the  interface  is  negligible.  However,  in  the  case  of  two 

rough  interfaces  the  correct  result  will  not  be  obtained  by  assuming 

the  two  interface  Rayleigh  reflection  coefficient  and  by  modifying 

4  <5 

it  to  include  the  local  incident  angles  as  has  been  done  by  Clay.  ’ 

The  only  other  paper  to  treat  scattering  by  rough  layers  is  by  Krishen 

and  Koepsel.^  They  recognize  the  inadequacy  of  Clay's  method  and 

attempt  to  account  for  the  scattering  between  layers.  They  assume 

7 

that  each  layer  is  in  the  farfield  of  the  other  and  use  Eckart's 
method  to  obtain  the  scattered  pressure  at  one  interface  due  to  the 
presence  of  the  other  rough  interface.  Further,  they  assume  that  the 
Rayleigh  reflection  coefficients  are  valid  at  each  interface  and  pro¬ 
ceed  on  the  basis  of  this  assumption.  Several  objections  can  be 
raised  to  their  method.  First,  Eckart's  method  does  not  lead  to  the 
correct  expression  for  the  pressure  scattered  or  transmitted  by  a 
rough  surface.  Further,  for  the  reasons  mentioned  previously,  it  is 
not  intuitively  obvious  that  the  Rayleigh  reflection  coefficients 
should  be  applied  to  the  scattering  at  each  interface. 

T Yr.  c..  assure  at  some  receiver  A  above  a  rough  surface 
is  given  by  the  potential  integral  in  Eq.  (4).  In  the  case  where  the 
water-bottom  interface  is  acoustically  penetrable  and  contains 
another  interface  at  some  deeper  depth,  the  pressure  density,  4,  at 
the  water- sediment  interface  must  include  the  interaction  and  inter¬ 
ference  effects  due  to  the  sub-bottom  reflections.  The  geometry  is 
shown  in  Dwg.  AS-71-354. 
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The  method  used  to  calculate  the  equivalent  pressure 
density  due  to  a  double  interface  is  identical  to  the  method  used  for 
the  single  interface,  but  in  addition  it  will  rely  heavily  on  the 
calculation  of  the  scattered  pressure  given  in  previous  final  reports. 
This  method  allows  the  integration  over  tne  lower  scattering  interface 
to  be  bypassed. 


Assume  that  the  two  interfaces  are  represented  by  the 
stochastic  processes  z^  =  C^(x,y)  and  z^  =  £g(x,y)  and.  that  the  mean 
planes  of  the  two  processes  are  parallel  and  separated  by  the  vertical 
distance  h.  Assume  further  that 

.MAX  .MAX  .  , 
h 

2 

and  that  the  insonified  area  of  the  surface  is  much  greater  than  h  and 
also  includes  many  correlation  lengths.  From  these  conditions  it  may 
then  be  assumed  that  the  pressure  density  (or  normal  velocity  field 
density)  at  any  point  on  the  lower  interface  is  the  sum  of  the  contribu¬ 
tions  from  a  region  of  the  upper  interface  which  is  many  correlation 
lengths  in  size.  The  significance  of  these  assumptions  is  that  now 
the  equivalent  pressure  density  must  be  calculated  in  terms  of  the 
average  pressure  densities  and  normal  velocity  densities  contributed 
by  the  lower  interface;  that  is,  the  averages  are  taken  before  calcu¬ 
lating  the  equivalent  surface  potential.  In  general,  taking  the 
averages  of  the  expressions  used  to  derive  the  equivalent  pressure 
density  is  much  easier  than  taking  the  average  of  the  equivalent  pres¬ 
sure  density  itself.  It  should  be  recognized,  however,  that  if  the 
contributing  surfaces  have  dimensions  of  only  fractions  of  a  correla¬ 
tion  length,  then  the  equivalent  pressure  density  must  be  calculated 
first  and  the  average  of  this  quantity  taken. 
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Using  these  assumptions,  the  equivalent  pressure  density 
in  the  case  of  a  double  interface  is  found  to  be 


where 


U  =  RPi  , 


R12  *  ay 


-^k^sinSgh 


1  +  ARi2R85e 


-i2kgSin02^ 


'ine  quantities  R.  „  and  R0„  are  the  reflection  coefficients  which  would 

Ld  dj 

exist  at  interfaces  1  and  2,  respectively,  if  medium  2  were  infinite  in 
extent;  they  are  given  by  Eq.  (17).  R2j  is  calculated  using  the  average 
pressures  at  interface  2.  Finally,  A  is  the  term  which  arises  because 
part  of  the  pressure  at  interface  1  has  been  scattered  upwards  by  the 
lower  interface. 


If  it  is  assumed  that  the  upper  interface  has  mild 
slopes,  then  A  is  given  by 


4- 


-ill  k  hcos0o+iiloC  ’k„cos0o> 
r)2cot  92)  e  ^  d  d  d  d 


where  tj2  is  the  local  slope  of  the  lower  interface,  and  is  the 
vaiiation  of  the  separation  of  the  interfaces  about  the  mean  value  h. 
Performing  the  average  indicated  in  the  previous  expression  yields 


2,  2  2*  V2 

s  k2cos  Sgh 


'  s  k?cos0  h  \ 

1  '  1  COt  @2  7  2  2  2  2  I?/2)  6 

v  |l  +  s  a  k2cos  02  / 


!(l+s2a2k2 


cos2@2) 


,  (20) 
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where  s  is  the  rms  value  of  the  lower  interface  slopes,  and  a  is  the 
rms  value  of  the  fluctuation  of  the  interface  separation  £ '  about  the 
mean  value  h. 


Clearly,  in  the  limit  of  two  plane,  parallel  interfaces 
A  ->  1  (since  s  and  a  -*  0)  and  R  approaches  the  Rayleigh  reflection 
coefficient  for  the  two  interface  case.  However,  it  must  be  noted 
that  Eq.  (19)  can  not  be  obtained  by  simply  taking  the  Rayleigh  reflec¬ 
tion  coefficient  and  referencing  it  with  respect  to  the  local  geometry. 
To  obtain  a  valid  expression  for  the  equivalent  pressure  density  in 
the  case  of  a  double  interface,  it  is  necessary  to  start  with  the 
appropriate  pressure  expressions  and  boundary  conditions. 

It  should  be  emphasized  that  R  [Eq.  (l9)]is  not  a 
constant,  but  rather  is  a  function  of  the  local  geometry  of  the  upper 
interface.  This  dependence  is  contained  in  the  R^2  terms,  which  are 
simply  the  reflection  coefficients  for  interface  1  referenced  to  the 
local  geometry.  Consequently,  R^  is  a  function  of  the  slopes  of  the 
upper  interface  and  must  be  used  in  conjunction  with  the  scattering 
integral  exactly  as  the  single  interface  reflection  coefficient  shown 
in  Eq.  (20). 


c.  Summary 

In  this  portion  [Section  B]  of  Chapter  II,  the  effect 
of  the  surface  configuration  and  physical  composition  on  the  boundary 
value  (i.e.,  the  equivalent  pressure  density)  used  in  the  scattering 
integrals  has  been  calculated.  This  work  was  motivated  by  the  need 
to  determine  the  boundary  values  in  a  manner  which  was  consistent  with 
the  integral  equation  of  scattering  and,  consequently,  required  no 
assumptions  about  the  direction  of  propagation  of  the  reflected  or 
scattered  field.;  For  the  case  of  a  single  interface,  it  was  found 
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that  if  the  source  was  many  wavelengths  from  the  surface,  the  correct 
boundary  value  for  a  rough,  acoustically  penetrable  surface  is  given 
by  multiplying  the  incident  pressure  density  by  the  Rayleigh  reflection 
coefficient  referenced  with  respect  to  the  local  slope.  In  the  event 
there  is  another  rough  interface  underlying  the  water-bottom  interface, 
the  equivalent  pressure  density  on  the  bottom  must  include  the  scat¬ 
tering  and  interference  effects  of  the  sub-bottom  interface.  However, 
it  was  found  that  these  effects  could  not  be  accounted  for  by  simply 
referencing  the  two  interface  Rayleigh  reflection  coefficient  to  the 
local  slope  because  this  does  not  account  for  the  scattering  by  the 
sub-bottom. 


The  scattering  and  interference  effects  of  the  sub-bottom 
can  be  properly  accounted  for  in  the  context  of  the  potential  theory 
while  simultaneously  avoiding  any  restrictions  on  the  direction  of 
propagation  of  the  waves  scattered  at  the  water-bottom  interface. 

The  single  and  double  interface  pressure  density 
coefficients  (reflection  coefficients)  given  by  Eqs.  (17)  and  (19) 
depend  upon  the  local  slopes  of  the  scattering  surface.  Hence,  these 
expressions  must  be  retained  in  the  scattering  integral,  and  averaged 
over  the  interface  slopes  along  with  the  other  slope  dependent  terms. 

The  remainder  of  Chapter  II  is  devoted  to  methods  of  calculating  the 
averages  over  the  surface  slopes  when  the  scattering  interface  is 
acoustically  penetrable.  In  addition,  since  the  slooes  play  an  impor¬ 
tant  role  in  both  the  scattering  and  the  acoustic  penetration  of  an 
interface,  considerable  attention  is  given  to  the  interaction  of  these 
two  effects  on  the  basis  of  their  slope  dependence. 

5*  Series  Evaluation  oi  the  Rayleigh  Coefficient 

In  order  to  treat  surfaces  with  fairly  large  rms  roughness, 
an  evaluation  of  the  reflection  coefficient  within  the  integral  for 
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the  scattered  pressure  was  carried  out  for  a  single  interface.  Using 
this  approach,  the  conditions  under  which  the  reflection  coefficient 
can  be  separated  from  the  scattered  pressure  integral  were  examined. 
The  series  expression  which  results  from  this  method  was  evaluated  for 
various  roughnesses  and  for  the  same  four  bottom  samples  used  in  the 
evaluation  of  Kuo's  reflection  coefficient. 

The  scattered  pressure  at  point  A  is  given  by  the  integral 


where 

M-  is  pressure  density  at  the  variable  point  N  on  the  surface, 
n  is  the  outward  normal  to  the  surface,  and 
r^  is  defined  in  Dwg.  AS-68-1100, 

The  value  of  4  on  the  boundary  is  given  by 

4(N)=~Rp.(N)  ,  (22) 


where  p^  is  the  value  of  the  incident  pressure  on  the  surface  and  R  is 
the  reflection  coefficient  given  by  Eq.  (17)  or  Eq.  (19),  which  will 
later  be  shown  to  be  a  function  of  the  local  slope  of  the  surface.  If 
it  is  now  assumed  that  the  incident  pressure  is  of  the  form 


ikr 


o 
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then  it  is  obvious  that 


H(N)  = 
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It  is  shown  in  the  Final  Report  under  Contract  N00024-69-C-1275  that 
this  integral  becomes 


PS(A) 


-//D^ 


e-ik7^  9^  9  -e  ) 
\  x  x  y  z ) 


e1  dxdy  ,  (26) 


where 


7  =  sinS,  +  sin0  , 
1  r 


and  is  a  unit  vector  directed  along  line  R^. 


If  it  is  assumed  that 


*  *  r\ 

e  •  e.  =  0 
X  1 


fe  .  e  =  cose 


'y  i 


e  •  e  =  sine 
z  1  r 


26 


then 


fA* 


A  A 

e  -e 

y  y  z 


)  •  e  =  -sin©  +  f  cos© 

/  1  r  Jy  r 


The  integral  for  the  scattered  pressure  becomes 


ps<A>  -*r}J  Do\— — 

S 


-iky? 


|^R(-sin©r+C 


ycos©r)Jix;dy  .  (30) 


At  this  point  a  suitable  expression  for  the  reflection 
coefficient  R  must  be  given.  Referring  to  Dwg.  AS-70-795>  p.  H,  the 
single  interface  reflection  coefficient  is  given  by 


p0c.smt,  -  p_csin<t0 

r(p  )  =  LA . 1 _ l  L-.  J. 

12  pgC^sinO^  +  p^c^sin<I>2 


■where  Eqs.  (15)  and  (16)  have  been  substituted  into  Eq.  (17)>  k  has  been 
replaced  by  and  =  0^+P  and  S>2  =  ^2  f  the  standard 

Rayleigh  reflection  coefficient  referenced  to  an  infinite  plane  oriented 
at  an  angle  0  to  the  mean  plane.  Application  of  Snell's  Law  to  the 
formula  is  accomplished  through  the  substitution 


sin$2  =  Vl  -  cos^&2  =  -  cos'’<J>^  , 


where  N  =  —  is  the  acoustic  index  of  refraction.  Applying  Eq.  (32) 
c2 


to  Eq.  (3l)  gives 


R(0l+P)  = 


pi  ci  rr 

p^sin^+p)  -  -i-i  \N 


N  -  cos  (0  +p) 


P,C  I  o  p 

P2c2sin(©L+p)  +  -  cos  (0^) 
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Expanding  sin(6^+p)  and  coc(0^+p)  and  using  various  trigonometric 
identities,  Eq.  (33)  ma;y  be  arranged  into  the  form 


R(e1,Ti)  - 


A 


A 


+  BT]  -  CyPT]2  +  Eq  +  F 
+  Bq  +  cVdt]1"  +  Eq  +  F 


(34) 


where 

7]  =  tanf3  ,  (35) 

A  =  PgCgSinO^  , 

B  =  poc,.cos0  , 

P1C1 

C  =  T1  •  (56) 

2  2 

D  =  N  -  sin  0^  , 

2 

E  =  sin  0 ,  and 

2  ? 

F  =  N  -  cos  0^ 


Since  q  is  the  tangent  of  P,  it  is  obvious  from  Dwg.  AS-70-795  thau  P 
is  just  the  slope  of  the  local  tangent  plane  with  respect  to  the  mean 
plane.  Because  the  formulation  is  restricted  to  the  plane  of  inci¬ 
dence  (ex  •  =  0).  it  may  be  assumed  that  q  =  Cy-  Using  this  result 

in  the  original  pressure  integral  gives 


(37) 

where  R(01,q)  is  given  by  Eq.  (34). 
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It  is  now  obvious  why  the  reflection  coefficient  cannot  be 

separated  from  the  integral  for  the  scattered  pressure  a  priori. 

The  dependence  of  R  upon  H,  a  quantity  which  appears  within  the  integral, 

places  restrictions  on  the  separation.  Here  it  is  recognized  that 

£  and  T]  are  random  variables  and  that  p  (a)  must  be  replaced  by  the 

s 

ensemble  average  <p  (A)>.  If  it  is  assumed  that  the  surface  height,  £, 

s 

and  the  surface  slope,  R,  are  independent,  then  the  integral  for 
<p  (A)>  takes  on  the  simple  form 


<P  (A)> 

s 


ik(R  +R  )' 
o  1 


- 01,n  )(  -  sinS^+Ticose^  dxdy 


(38) 


To  separate  an  effective  reflection  coefficient  (in  an  averaged  form) 
from  the  integral  it  must  also  be  assumed  that  T)  is  stationary  so  that 


Since  f)  is  also  stationary,  its  expectation  value  may  also  be  removed 
from  the  integral,  leaving  only  a  form  which  is  immediately  recognized 
as  the  pressure  reflected  from  a  plane  surface,  p^, 


<Ps(A)>=  -[<R(ei,T1)>  -  coter<T1R(ei,T])>]  <eliQ^>pr 


(l+o) 
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Assuming  a  Gaussian  distribution  of  heights,  the  average  over  £  gives 


<Ps(A)>  =  -  [<R(©1,il)>  -  eot0r<TiR(01,ri)>J  e_g/2pr  .  (4l) 

2 

where  g=(k7cr)  and  c  is  the  rms  height  value. 

The  expression  in  brackets  in  Eq.  (4l)  may  be  regarded  as 
an  effective  reflection  coefficient  which  has  been  averaged  along  with 
the  slope  terms  in  the  integral  for  the  scattered  pressure  and  then 
separated  from  the  integral  in  an  averaged  form.  The  expectation  values 
<R>  and  <t]R>  are  still  quite  difficult  to  evaluate  directly  due  to  the 
complex  form  of  R.  In  order  to  evaluate  these  expressions,  R(  9^,t)  ) 
is  expanded  in  a  Maclaurin  series; 


^0 


The  first  term  of  this  series  is  equivalent  to  the  Rayleigh  reflection 
coefficient  for  a  plane  surface.  This  series  expansion  of  R(0^,q)  has 
been  widely  used  in  rough  surface  scattering  literature  (Refs.  8-lU). 
Using  this  expansion  the  evaluation  of  the  reflection  coefficient  is 
straightforward  and  is  given  by 

|<R>  -  coter<T)R>J  =  £r(1  -  <T}>  cote^)  +  R'  (<T)>  -  <r\^>  cot ©r) 

+  —  (<rj‘>  -  <rj^>  cote^)  +  .  .  .J  >  (^) 


where  R  and  its  derivatives  are  evaluated  at  q=0. 


The  only  expectation  values  to  be  calculated  are  of  the  form 


<T)  >  = 


co 

/ 


T)  P(lj)dT)  , 


(^5) 


-tan©. 


where  the  lower  limit  of  -tan©^  represents  a  simple  shadowing  considera¬ 
tion  for  backscattering  geometry.  In  forward  and  specular  scattering 
the  upper  limit  is  replaced  by  the  tangent  of  the  receiver  angle. 

For  a  Gaussian  distribution  of  surface  slopes  and  for  l,  an  odd  integer, 
Eq.  (45)  becomes 
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where  s  is  the  rms  value  of  the  surface  slope  and  P  is  the  incomplete 
gamma  function.  Similarly  for  l  even, 


For  the  forward  and  specular  scattering  cases  the  expectation  values 
of  h  contain  additional  incomplete  gamma  functions  depending  on  6  . 

These  expectation  values  may  be  easily  calculated  and  combined  with 
the  values  of  R,  R',  R",  etc.,  to  find  the  value  of  the  effective 
reflection  coefficient  [<R>  -  cot6?r  <RR>1 . 

It  should  be  noted  here  that  the  expression  [<R>  -  cotf^  <t)R>], 
which  has  been  called  the  effective  reflection  coefficient,  actually 
contains  terms  which  arise  from  the  exact  treatment  of  the  surface 
slopes  in  the  scattered  pressure  integral  as  well  as  from  the  penetra¬ 
bility  of  the  surface.  The  rigorous  treatment  of  the  surface  slopes 
produces  the  term  (-sin@r+  qcose^)  in  the  integrand  of  Eq.  (30).  This 
term  is  a  result  of  the  scattering  process  and  is  independent  of  the 
penetrability  of  the  surface.  Averaging  this  term  together  with  R  is 
necessary  because  they  both  depend  on  q,  but  this  averaging  results  in 
an  effective  reflection  cofficient  in  which  effects  of  scattering  are 
combined  with  effects  of  penetrability.-  To  see  this  more  clearly,  it 
is  only  necessary  to  examine  the  behavior  of  the  effective  reflection 
coefficient  for  an  impenetrable  rough  surface.  For  a  pressure  release 
surface  R(©^,q)  =  -1,  so  the  effective  reflection  coefficient  is  given 

by 


[l  -  cot&  <rj  >] 
r 
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If  <t]>  is  calculated  according  to  Eq.  (4l)  with  the  appropriate  limits 
on  the  integral  for  each  scattering  geometry,  tie  results  are  given  t>y 


[l  -  COt0  <Tj>] 
r 


SPECULAR 


=  1 


J 


[l  -  c °t  0r<Ti >]  BAC  Kgc ATTER 


tan  0  1 


and 


(IS) 


[l  -  cot0r<ri>] 


FORWARD  SCATTER 


These  results  are  shown  for  various  values  of  S(rms  slope)  in 
Dwg.  AS-7I-294  for  backscattering  geometry  and  in  Dwg.  AS-71-295  for 
forward  scattering  geometry.  This  expression  represents  the  contribu¬ 
tion  to  the  effective  reflection  coefficient  that  is  due  solely  to  the 
scattering  process.  In  Dwg,  AS-71-295>  curves  corresponding  to  the 
larger  values  of  s  were  not  plotted  because  they  all  lie  quite  close 
to  the  s  =  0.5  curve.  Realistic  values  of  the  rms  slope  of  the  ocean 
bottoms  are  generally  quite  small.  For  the  model  surfaces  used  in  this 
program  the  rms  slopes  range  from  s  =  0.1  to  s  =  0.3.  From 
Dwgs.  AS-71-294  and  AS-71-295  it  is  obvious  that  for  these  small  slope 
values,  the  slope  dependent  term  is  very  nearly  unity  except  at  very 
low  grazing  angles.  Since  this  scattering  contribution  can  be  calcu¬ 
lated  exactly,  a  true  reflection  coefficient,  which  depends  only  upon 
the  penetrability  of  the  bottom,  may  be  defined  as 


_  [<R>  -  cot0^<qR>] 

R  =  . ll  “  COt0r<q>j 


(49) 
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SLOPE  DEPENDENT  TERM  (0  -  cot  0r  <^>])  OF  THE  SCATTERED 
PRESSURE  FOR  AN  IMPENETRABLE  SURFACE  AND  BACKSCATTERING 
GEOMETRY  FOR  VARIOUS  VALUES  OF  THE  rms  SLOPE 


■where  the  expression  in  the  numerator  is  the  effective  reflection 
coefficient  and  the  expression  in  the  denominator  has  been  shown  to 
be  the  scattering  contribution  to  the  effective  reflection  coefficient. 
Using  this  result  the  scattered  pressure  for  a  Gaussian  distribution  of 
heights  if  given  by 

p  =  R  [l  -  cot0  <R>]  e“g/2  p  , 
s  r  r 

where  R  depends  only  upon  the  penetrability  and  obviously  is  unity 
for  an  impenetrable  or  pressure  release  surface.  Because  the  slope 
dependent  term  [l  -  cot0r<rj>]  is  nearly  unity,  for  the  cases  considered 
here  it  will  be  assumed  that  the  effective  reflection  coefficient  pro¬ 
vides  an  adequate  approximation  to  the  true  reflection  coefficient 
given  by  Eq.  (1+9). 

The  effective  reflection  coefficient  [<R>  -  cot0^<r)R>] 
will  now  be  evaluated  for  a  penetrable  surface  under  the  assumption 
that  it  adequately  approximates  the  true  reflection  coefficient  and 
thus  represents  the  change  in  the  scattered  pressure  due  to  the  pene¬ 
trability  of  the  bottom. 

To  allow  for  attenuation  in  the  bottom  the  velocity  c^  must 
be  complex  with  an  imaginary  part  that  is  proportional  to  the  attenua¬ 
tion.  This  makes  the  constants  A,  B,  C,  D,  and  F  in  the  expressions 
for  R,  R',  R",  etc.,  complex,  and  hence  the  entire  effective  reflection 
coefficient  is  a  complex  quantity.  Plots  of  the  magnitude  of  this 
coefficient  for  various  bottoms  with  several  values  of  the  rms  slope 
are  given  in  this  section. 

Drawings  AS-71-24,  AS-71-25,  and  AS-71-26  show  the  reflection 
coefficients  for  various  source  and  receiver  geometries  for  Mackenzie's 
sample  number  J2  at  a  frequency  of  5  kHz.  The  measured  attenuation 
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REFLECTION  COEFFICIENT  vs  INCIDENT  GRAZING  ANGLE 
FOR  SPECULAR  SCATTERING  GEOMETRY 
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of  4.48  dB  per  ft  at  3°  kHz  wee  extrapolated  to  kHz  assuming  that 
attenuation  varied  as  the  square  root  of  the  frequency.  Drawing  AS-71-24 
shows  the  magnitude  of  the  effective  reflection  coefficient  plotted 
versus  grazing  angle  for  a  specular  scattering  geometry  (0  =  8  ). 

In  the  limit  of  no  roughness  the  first  term  in  the  series  expansion 
of  R(0,q)  is  the  only  nonzero  term.  This  term  is  just  Mackenzie's 
original  result  for  a  plane  attenuating  surface  and  is  plotted  as  a 
solid  line  in  Dwg.  AS-71-24.  The  dashed  line  represents  a  small  rms 
slope  of  0.05  and  shows  little  deviation  from  the  zero  slope  term. 

The  dotted  line  represents  a  larger  rms  slope  value  of  0.10  and  exhibits 
large  fluctuations  in  the  region  around  the  critical  angle  which  occurs 
at  a  grazing  angle  of  30  deg  for  this  surface.  The  failure  near  the 
critical  angle  is  a  problem  inherent  in  the  formulation  of  this  method 
and  will  he  discussed  later.  Similar  results  are  obtained  for  back- 
scatter  and  are  shown  in  Dw^,.  AS-71-25-  Here  the  plane  surface  term 
has  no  meaning  and  is  not  plotted.  The  forward  scatter  curves  shown 
in  Dwg.  AS -71-26  are  typical  and  have  no  fluctuations  as  long  as  the 
incident  grazing  angle  is  not  in  the  region  immediately  around  the 
critical  angle.  Similar  results  are  shown  in  Dwgs.  AS -71-214  and 
AS-71-215  for  bottom  sample  number  28,  which  exhibits  a  critical  angle 
near  If  deg.  Hie  forward  scatter  curves  were  not  plotted  since,  for 
these  rms  roughness  values,  the  forward  scatter  coefficient  was  very 
nearly  a  constant  (R=0.29).  Samples  number  ?7  and  13  do  not  have  a 
critical  angle  but  instead  have  an  angle  of  intromission  at  about 
6  deg.  The  reflection  coefficient  curves,  shown  in  Dwgs.  AS-71-216, 
AS-71-217,  AS-71-218,  and  AS-71-219,  are  relatively  smooth  because  the 
angle  of  intromission  is  actually  below  the  lower  limit  of  the  range  of 
grazing  angles  of  interest.  It  is  expected  that  an  angle  of  intro¬ 
mission  will  produce  the  same  kind  of  fluctuacions  in  the  reflection 
coefficient  as  that  observed  near  the  critical  angles ,  Again  the 
forward  scatter  curves  were  omitted  because  they  were  very  nearly 
constant. 
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It  was  first  thought  that  the  fluctuations  near  the  critical 

angle  were  simply  the  result  of  not  retaining  an  adequate  number  of 

terms  in  the  series  expansion  given  by  Eq.  (25).  However,  a  partial 

sum  decomposition  of  the  series  shows  that  this  is  not  the  case  and 

that  the  series  is  actually  diverging  in  this  region.  This  indicates 

that  the  problem  is  inherent  in  the  formulation  of  the  method.  The 

presence  of  fluctuations  near  the  critical  angle  points  out  several 

similarities  between  this  approach  and  the  method  used  by  Brekhovskikh 

to  treat  the  reflection  of  spherical  waves  from  a  plane  interface, 

15 

which  suffers  similar  difficulties.  Brekhovskikh  uses  a  plane  wave 
decomposition  of  a  spherical  incident  wave,  and  saddle  point  integra¬ 
tion  to  arrive  at  an  effective  reflection  coefficient  given  by 


where 

H  =  |  [r"(90>  +  B'(90)c°teoJ  . 

Here  R'  (9  )  and  R"($  )  are  the  derivatives  of  the  reflection  coef- 
o  o 

ficient  with  respect  to  9.  The  reflection  coefficient  itself  is  of 
the  same  form  as  Eq.  (17).  In  the  derivation  of  this  spherical  wave 
reflection  coefficient,  it  is  assumed  that  R(0)  is  slowly  varying. 
Near  the  critical  angle  this  is  no  longer  true  and  the  derivatives  of 
R  with  respect  cc  9  become  large.  Thus,  as  Brekhovskikh  shows,  the 
series  is  no  longer  valid  near  the  critical  angle. 

The  application  of  these  ideas  to  the  series  given  by 
Eq.  (27)  is  straightforward.  It  was  initially  assumed  that  the  first 
term  of  Brekhovskikh ' s  series  was  an  adequate  approximation  to  the 
spherical  wave  reflection  coefficient.  This  approximation  is  quite 
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good  because  kR^  is  very  large  for  the  scattering  geometries  being 
considered.  At  this  point  the  failure  of  the  series  development 
presented  in  this  section  is  difficult  to  understand  because  the  terms 
that  cause  Brekhovskikh ' s  coefficient  to  be  invalid  near  the  critical 
angle  have  simply  been  neglected  on  geometrical  grounds.  However,  when 
the  effects  of  roughness  are  included,  the  reflection  coefficient 
takes  on  the  form  R(@+0),  where  3  is  the  random  angle  associated  with 
the  slope  of  the  local  tangent.  In  order  to  evaluate  R  it  is  necessary 
to  expand  it  in  the  Maclaurin  series  given  by  Eq.  (25).  Since  the 
reflection  coefficient  is  of  the  form  R(0+p),  it  can  easily  be  shown 
that 


dR  dR  dR  dn  dR  20 

Z9  =  =  3n  §  =  sec  P 

Thus,  the  derivatives  of  R  with  respect  to  q  which  appear  in  the 
Maclaurin  series  behave  in  much  the  same  manner  as  the  derivatives  of 
R  with  respect  to  9  which  had  previously  been  eliminated  from 
Brekhovskikh ' s  formula.  Consequently,  the  effective  reflection  coef¬ 
ficient  given  by  Eq.  (27)  suffers  the  same  fluctuations  at  grazing 
angles  near  the  critical  angle.  If  Eq.  (24)  could  be  evaluated  directly, 
then  the  fluctuations  associated  with  the  expansion  could  be  avoided. 
Direct  evaluation  of  Eq.  (2k)  might  be  accomplished  using  numerical 
integration  techniques  although  this  approach  is  complicated  by  the 
fact  that  R( 0  ,q)  is  a  complex  function  of  real  arguments.  An  alter¬ 
nate  method  utilizing  the  complex  nature  of  R  would  be  to  express 
<2R( 0^ , iq )>  and  <qP(0^,q)>  as  contour  integrals.  The  evaluation  of 
these  contour  integrals,  however,  is  not  straightforward  because  of 
the  presence  of  branch  cuts.  These  ideas  are  being  investigated  at 
the  present  time. 

To  summarize  the  results  of  this  subsection,  it  was  shown 
that  an  effective  reflection  coefficient,  for  a  single  interface,  which 
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could  be  removed  from  the  integral  for  the  scattered  pressure  in  an 
averaged  form  -was  given  by 

|<R>  -  coter 

It  was  also  shown  that,  for  small  values  of  the  rms  slope,  the 
scattering  contribution  to  this  expression  could  be  neglected  and  that 
the  effective  reflection  coefficient  represented  the  change  in  the 
scattered  pressure  due  to  the  penetrability  of  the  surface.  For 
larger  slope  values  the  true  reflection  coefficient  was  given  by 
Eq.  (49)  and  can  be  easily  calculated  from  the  effective  reflection 
coefficient.  In  order  to  calculate  the  effective  reflection  coeffi¬ 
cient  for  penetrable  rough  surfaces,  R(0^,q)  was  expanded  in  a 
Maclaurin  series  in  tj.  This  technique  made  the  evaluations  of  <R> 
and  <t]R>  quite  simple  but  introduced  fluctuations  near  the  critical 
angle,  or  angle  of  intromission.  It  was  shown  that  these  fluctuations 
were  a  result  of  using  the  Maclaurin  expansion  in  a  region  where  the 
function  was  not  slowly  varying.  The  closed  form  expression  for  the 
effective  reflection  coefficient,  [<R>  -  cot0^  <tjR>],  is,  however,  the 
result  of  a  rigorous  derivation,  and  so  the  Maclaurin  expansion  can 
be  expected  to  give  the  correct  results  in  the  regions  where  R  is  slowly 
varying.  This  gives  a  reflection  coefficient  which  is  valid  over  all 
angles  for  surfaces  where  R  varies  slowly  and  is  valid  over  a  large 
range  of  angles  for  surfaces  that  exhibit  a  critical  angle  or  an  angle 
of  intromission.  In  addition,  if  the  rms  slope  is  very  small,  as  it 
often  is  for  realistic  ocean  bottoms,  the  fluctuations  in  the  derivatives 
of  R  in  the  Maclaurin  series  are  insignificant  because  the  derivatives 
are  multiplied  by  expectation  values  of  powers  of  tj  which  are  extremely 
small.  Thus,  good  results  can  be  obtained  near  the  critical  angle  if 
the  slopes  are  small  enough.  This  aspect  is  illustrated  by  Dwg.  AS-71-24, 
where  good  results  are  obtained  for  s=0.05  but  fluctuations  occur  when 
s=0 . 1 . 
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The  results  given  in  this  subsection  are  derived  specifically 

for  a  single  interface.  The  double  interface  problem,  however,  can  be 
treated  in  a  similar  manner  and  will  also  result  in  an  expression  that 
fluctuates  near  the  critical  angle  because  the  expression  will  still 
involve  the  derivatives  of  R^.  It  is  expected  that  the  effective 
reflection  coefficient  for  a  double  interface  will  be  valid  over  approxi¬ 
mately  the  same  range  of  angles  and  slopes  as  the  single  interface 
reflection  coefficient.  However,  it  should  again  be  pointed  out  that 
these  expressions  are  valid  for  a  large  number  of  the  combinations  of 
slopes,  angles,  and  bottom  parameters  of  practical  interest. 

4.  Stationary  Phase  Evaluation  of  the  Reflection  Coefficient 

As  seen  in  the  previous  subsection  on  the  series  evaluation 
of  the  refie.  tion  coefficient  for  the  single  interface  case,  accurate 
re:-’.! '  ts  can  be  obtained  for  certain  special  cases.  In  the  event  that 
•he  source  crazing  angle  is  near  the  critical  angle  (or  angle  of  intro¬ 
mission,  depending  on  whether  the  index  of  refraction  N  is  greater  than 
or  les  than  l)  and  chat  the  rms  slope  is  moderate,  the  series  evalua¬ 
tion  fai.'.c  because  of  the  behavior  of  the  derivatives  of  R.  A  series 
evaluation  of  the  double  interface  reflection  coefficient  would  fail 
for  identical  reasons. 

A  method  of  evaluating  the  reflection  coefficient  which  is 
valid  for  the  case  where  both  the  rms  heights  and  slopes  are  moderate 
or  large  is  given  by  the  stationary  phase  evaluation.  This  method 
avoids  the  difficulties  encountered  in  the  series  method.  In  the  sta¬ 
tionary  phase  method  the  reflection  coefficient  is  evaluated  at  the 
angle  0  which  is  the  local  grazing  angle  at  the  stationary  phase  points. 
If  <p  is  the  angle  between  the  plane  of  incidence  and  the  scattering 
plane,  the  assumption  cp=0  implies  that  only  the  scattering  in  the 
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plane  of  incidence  will  be  of  interest.  The  local  angle  0  is  then 
given  by 


sin'J  - 


sin0.  -  £  cose. 
1  x  1 

/  2  2 

(1+5x+V 


(51) 


where  at  the  stationary  phase  points 


v  v 

C  =  _JL  {-  , 

v  v  ’  x  v 

v  7  rr 


and 


v^  =  k(cos0^  -  cosQ^cosep)  =  kCcose^ 

v^  =  -kfcose^sincp)  =  0  , 

v  =  -k(sine  +  sin0.) 
z  r  l 


cose  ) 

r 


This  gives 


1  l/° 

sin®  _  —  (l  -  cose,  cose  +  sine,  sine  )  /  -  , 

/2  1  r  1  r 


(52) 


/e.  +  e  \ 

=  sin("--r~~  |  , 


where  6^  is  the  source  grazing  angle  and  0^  is  the  receiver  angle.  If 
e^+p  instead  of  $  is  used  in  the  reflection  coefficients  given  by 
Eqs.  (17)  and  (19),  then  for  the  single  interface  case  the  reflection 
coefficient  becomes 


8(0^11)  -» R 


(53) 
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and  for  the  double  interface  case  the  reflection  coefficient  becomes 


(8.  +  8 
1  r 


R^,".) 


12\  2 


-i2k2sin0oh 


1  +  AR_R 


^i+  0  \  -i2k2sin9gh 


23  12V  2 


(5*0 


It  should  be  noted  that  R  given  by  Eq.  (53)  is  identical  to  the  R^2 
contained  in  Eq.  (5*0*  The  other  terms  in  Eq.  (5*0  are  constants 
with  respect  to  the  stationary  phase  calculation;  consequently,  in 
the  following  chapter  only  Eq.  (53)  "will  "be  discussed. 


For  the  backscatter  case,  0  =  it  -  0. ,  so  =  it/2.  Thus  at 

all  incident  grazing  angles,  the  reflection  coefficient  for  back¬ 
scatter  geometry  is  a  constant.  In  the  specular  direction,  0^  =  0^,  so 
$  =  0^.  This  simply  reproduces  Mackenzie’s  original  curves  for  the 
reflection  coefficient  from  a  smooth  penetrable  surface.  For  the 
forward  scatter  case,  <I>  =  (0^+0^)/2.  The  curves  for  this  case  are 
given  in  Dwg.  AS-71-28.  These  results  are  particularly  easy  to  inter¬ 
pret  in  that  each  point  on  the  surface  is  assumed  to  reflect  in  the 
specular  direction  with  respect  to  the  local  tangent  plane  at  that 
point.  Thus,  the  reflection  coefficient  depends  only  upon  the  scat¬ 
tering  geometry.  Once  the  position  of  the  receiver  is  specified,  the 
only  points  which  can  contribute  to  the  field  at  the  receiver  are  the 
stationary  phase  points. 

Examination  of  Eq.  (37)  reveals  that  the  results  given  in 
Eqs.  (53)  and  (5*0  are  not  yet  complete.  According  to  Eq.  (37),  the 
slope  dependence  is  contained  in  both  R  and  the  scattering  integrand  as 


R(9^,q)[sin©r  -  TjcosS^] 


(55) 
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Eq.  (55), 


Considering  the  complete  slope  dependence  contained  in 
the  stationary  phase  result  then  becomes 


/e.  +  e  \ 

R(0.  ,I))[sin6  -  T|cos0  ]  =  Rl-“ — -Jf(0.,0  ) 


1  r 


(56) 
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where  R  is  given  by  either  Eq.  (53)  or  Eq*  (54),  f-ed  F  is  Beckmann's 
F  and  is  given  bj 


1  -  cos (0.  d  ) 

F(9i’V  -  -sina:  -rkr-  ■ 

1  r 

Previously,  it  was  pointed  out  chat  because  both  R  and  the 
scattering  term  given  in  the  brackets  in  Eq.  (55)  are  slope  dependent, 
the  acoustic  penetrability  of  the  surface  interacted  with  the  scatter¬ 
ing.  This  same  effect  occurs  when  a  stationary  phase  calculation  is 
made  since  the  term  which  is  then  obtained  is  not  just  R  but  rather  RF 
given  by  Eq.  (56). 

The  stationary  phase  calculation  presented  here  is  a 

l6  IT 

generalization  of  the  results  of  Hagfors  and  Stogryn.  On  the 
basis  of  a  stationary  phase  calculation,  they  obtained  a  reflection 
coefficient  which,  in  the  case  of  backscattering,  was  simply  a  constant 
factor  equal  to  the  Rayleigh  reflection  coefficient  evaluated  at 
normal  incidence. 
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III.  THE  EFFECT  OF  A  ROUGH,  PENETRABLE  LIQUID  BOTTOM 
ON  SONAR  OPERATION 

The  motivation  for  the  theoretical  work  presented  in  the  previous 
section,  and  in  previous  final  reports,  has  been  to  develop  realistic 
models  of  sonar  operation  in  the  ocean  environment.  Prior  studies 
have  verified  the  mathematical  models  which  apply  to  an  ocean 
boundary  composed  of  plane  layered  sediments  or  a  pressure  release 
sinusoidal  surface.  More  recent  work  has  resulted  in  the  develop¬ 
ment  of  a  model  of  forward  and  backward  scattering  at  the  sea 
surface.  The  previous  theoretical  work  was  somewhat  simpler  than  the 
present  case,  since  the  scattering  effects  were  due  only  to  topog¬ 
raphy.  In  developing  a  model  for  a  rough  ocean  bottom,  an  additional 
complication  is  introduced,  since  now  both  the  topography  and  the 
physical  composition  of  the  bottom  affect  the  scattering  behavior. 

The  model  that  will  be  presented  in  this  report  is  incomplete  in 
several  aspects.  First,  the  model  is  restricted  to  liquid  bottom 
types;  that  is,  the  shear  mode  of  propagation  must  be  insignificant 
in  comparison  to  the  compressional  or  longitudinal  mode  of  propaga¬ 
tion.  Second,  this  model  has  not  as  yet  been  experimentally  verified. 

Special  aspects  of  the  simplest  case  (single  interface  case)  have 

17 

been  discussed  in  the  literature  (Hagforsx  and  Stogryn,  for 
example),  but  a  ge;  ->ral  verification  will  not  be  completed  until 
the  end  of  the  next  contract  year.  The  preliminary  experimental  work 
is  reported  in  Appendix  A,  however. 

The  effects  of  a  rough  bottom  are  encountered  in  a  number  of 
situations .  In  general  sonar  operation  the  presence  of  a  rough 
penetrable  bottom  will  degrade  the  effective  range  of  the  sonar  and 
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also  give  rise  to  a  reverberation  return  which  may  mask  the  presence 
of  a  target.  For  a  sonar  which  operates  in  a  bottom  bounce  mode  the 
type  of  bottom  present  in  the  area  of  operation  is  of  critical 
importance.  If  the  bottom  is  smooth  but  penetrable,  then  the  sonar 
system  will  suffer  only  a  loss  in  signal  level.  Where  the  bottom  is 
rough  but  relatively  impenetrable,  the  sonar  system  will  be  degraded 
in  several  respects.  The  scattering  effect  of  the  bottom  will  be 
observed  as  a  general  beam  broadening  with  consequent  loss  of  bearing 
accuracy.  Further,  a  rough  bottom  will  give  rise  to  fluctuations  in 
the  returned  signal.  Obviously,  a  rough,  penetrable  bottom  will  give 
rise  to  a  combination  of  the  effects  just  mentioned. 

It  is  important  to  determine  the  dominant  effect  of  the  bottom; 
that  is,  whether  penetration  or  roughness  is  the  dominant  loss  factor. 
The  importance  of  this  point  is  test  illustrated  by  an  example. 

Assume  that  the  perpendicularly  measured  bottom  loss  in  two  regions 
was  -20  dB.  If,  in  one  region,  this  loss  was  primarily  due  to 
roughness,  then  at  the  grazing  angles  for  which  a  bottom  bounce  sonar 
is  operated  the  loss  will  still  be  large  and  significant  signal  fluctua¬ 
tion  and  beam  degradation  will  also  occur.  Suppose  that  in  the  other 
region  the  loss  was  due  mainly  to  penetration.  Then  at  the  usual 
operating  angles  the  loss  will  perhaps  be  only  -5  dB  (assuming 
cB0TT0M  >  S/ATER^’  an<*  associ-a^e^  beam  broadening  and  signal  fluctua¬ 
tions  will  be  insignificant.  Consequently,  in  the  second  region  the 
bottom  bounce  mode  may  be  effectively  employed  at  low  grazing  angles. 
These  considerations  take  on  added  weight  when  it  is  realized  that  it 
is  impossible  to  acoustically  map  the  entire  ocean;  hence,  the  local 
sonar  operator  will  occasionally  have  to  make  a  determination  of  the 
suitability  of  the  bottom  for  the  bounce  mode. 


A.  Physical  Description  of  Some  Ocean  Bottoms 


To  be  worthwhile,  the  effects  of  the  bottom  presented  in  this 
report  should  represent  some  fairly  typical  bottom  types.  However, 
those  bottoms  which  are  composed  of  rock  or  compacted  sediments  must 
be  excluded  since  the  shear  mode  will  be  an  important  propagation 
mechanism.  Essentially,  these  regions  are  centered  near  the 
midoceanic  ridges  which  are  generally  considered  to  be  topographi¬ 
cally  unsuitable  for  bottom  bounce  operation  anyway. 

The  other  major  physiographic  provinces  that  are  of  interest, 
which  are  also  thought  to  be  usually  good  reflectors,  are  the  conti- 
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nental  margins  and  the  ocean  basin  floor.  The  continental  margins 
are  usually  found  at  moderate  depths  (100  -  1000  fathoms),  have  thick 
sediment  layers,  and  the  topography  may  consist  of  a  small  average 

1  ft 

gradient  (l:200  -  1:600)  with  small  local  relief  (10  -  20  meters). 

19 

An  even  smaller  scale  relief  may  be  due  to  strewn  rocks  and  boulders 

or  to  the  characteristic  sinusoidal  ridges  resulting  from  turbidity 

currents."*"  The  principal  reflections  from  the  continental  margin 

may  then  either  be  due  to  the  bottom  itself  or  to  a  sub-bottom  that  is 

21 

2  to  10  meters  deep.  These  features  are  characteristic  of  the 
Atlantic  Ocean  and  particularly  the  North  Atlantic  Ocean. 

The  ocean  basin  floors  are  generally  f_und  at  deeper  depths 

(1000  -  4000  fathoms),  have  a  nonexistent  or  thin  (0.2  -  4  mener) 

19 

organic  sediment  layer  overlying  a  fairly  dense  bottom,  and  the 

topography  may  consist  of  an  average  gradient  less  than  1:1000  and, 

except  for  isolated  hills  and  sea  mounts,  relative  relief  of  less 

18 

than  1  meter.  The  ocean  basin  is  thought  to  be  quite  smooth,  and 
where  it  is  not,  the  organic  sediment  tends  to  smooth  over  the  small 
scale  features.  The  principal  reflections  from  the  e ui  basin  floors 
are  due  either  to  the  bottom  or  a  very  near  lying  sub-bottom.  These 
features  are  characteristic  of  large  areas  of  the  Pacific  and 
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Indian  Oceana,  and  to  a  lesser  extent  to  certain  portions  of  the 
Atlantic  Ocean. 

In  addition  to  the  fact  that  the  physical  description  of  the 
ocean  bottom  is  well  correlated  with  the  type  of  physiographic 
province,  the  acoustic  parameters  are  also  moderately  well  correlated 
with  the  physiographic  provinces .  A  consequence  of  these  relation¬ 
ships  is  that  acoustic  domains  in  the  ocean  largely  coincide  with 
physiological  provinces. 

Table  I  gives  the  parameters  and  the  physical  description  of  the 

four  ocean  bottom  models  which  will  be  used  in  the  study  of  the  effect 

of  a  rough,  penetrable  bottom  on  bottom  bounce  sonar  operation.  The 

velocity  and  density  parameters  used  in  these  models  are  based  pri- 

3  22 

marily  on  the  data  presented  in  Hamilton,  et  al.,  Sutton,  et  al  , 
Shumway,  D  Hamilton,  Nafe  and  Drake,  and  Hampton  for  different 
physiographic  provinces  in  both  the  Atlantic  and  Pacific  OceanB. 

Note  that  in  the  two  layer  models  given  here,  the  upper  layer  is  a 
low  speed  layer;  i.e.,  c^pj,^  <  C-^ATER  •  The  shear  velocity  has  been 
omitted  since  it  is  assumed  that  the  bottoms  behave  as  liquid  layers. 

The  absorption  used  in  Table  I  is  based  on  the  data  in 

2  5  2*7  28  28 

Sbnmiway,  Wood  and  Weston,  Cole,  and  Hampton.  The  absorption 

is  based  on  an  assumed  first- power  frequency  dependence,  which  allows 

the  wave  number,  k,  to  be  wriiLen  as 


> 


where  the  absorption  per  wavelength,  a',  given  in  Table  I  is  then 
related  to  the  usual  absorption  a  given  in  the  sonar  equation  by 


a'  =  xa 
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The  advantage  of  writing  a'  instead  of  a  is  that  ex'  is  a  nondimens ional 
quantity,  and  if  the  frequency  does  obey  a  first- power  law,  then  the 
absorption  at  any  frequency  is  easily  obtained  using  the  simple 
relation  given. 

B.  The  Theoretical  Prediction  of  the  Forward  Reflection  and 
Scattering  of  a  Sonar  Beam  by  a  Penetrable  Bottom 

1  h-  5 

It  has  been  postulated  by  Kuo  and  assumed  by  Clay  *  that  the 
scattering  behavior  of  a  penetrable,  rough  bottom  is  given  very 
simply  by 


ofn  =  |R(ei)|2as  ,  (57) 

where  a  is  just  the  scattering  coefficient  for  ar.  impenetrable  rough 
s 

surface  and  is  defined  as  the  ratio  of  the  intensity  scattered  by  the 
impenetrable  rough  surface  in  the  direction  of  interest  to  the 
intensity  reflected  in  the  specular  direction  by  an  impenetrable 
plane  surface.  According  to  Clay,  R( )  is  the  Rayleigh  reflection 
coefficient  defined  at  the  given  incident  angle,  9 ,  for  perfectly 
plane  layers  with  parameters  identical  to  the  rough  layers.  Hence, 
R(@1)  is  defined  by  Eq.  (35)  with  the  local  slope  angle  0  set  to 
zero.  According  to  Kuo,  R(9  )  is  given  by  his  8(7)  [Eq.  (lj]  for  the 
case  of  backscattering. 

The  viewpoint  implicit  in  the  expressions  given  by  both  Kuo  and 
Clay  is  that  ihe  scattering  features  due  to  the  physical  px-operties 
of  the  bottom  and  the  scattering  features  due  to  the  acoustical 
properties  of  the  bottom  are  completely  separable  and  may  be  inde¬ 
pendently  expressed  by  the  terms  c  and  R( 0  ),  respectively.  In 
general  this  viewpoint  is  incorrect  since  both  R(@^)  and  the  scatter¬ 
ing  factor  depend  on  the  local  rough  surface  slope  value .  These 
terms  must  be  combined  when  the  average  over  the  slopes  is  to  be  taken. 
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When  extended  to  multiple  layer  media,  the  breakdown  of  Clay's 
viewpoint  is  even  more  severe  since,  as  seen  in  Eq.  (19),  the  reflec¬ 
tion  coefficient  now  depends  on  both  the  local  slope  and  the  local 
separation  of  the  layers  [which  depends  on  the  local  height  £(x,y)]. 
Hence,  the  effects  of  the  physical  and  acoustical  properties  of  the 
bottom  are  now  bound  together  in  both  the  average  over  the  surface 
slopes  and  the  average  over  the  surface  heights. 

These  arguments  may  be  mitigated  under  several  circumstances, 

but  only  for  the  case  of  a  single  interface.  If  the  slopes  are 

small,  then  for  many  scattering  configurations  they  may  be  ignored. 

Or  in  cases  where  the  slopes  are  significant,  they  may  still  be 

ignored  in  forward  scattering  if  the  grazing  angles  are  noc  too 

small.  In  these  two  cases  Eq.  (57)  as  given  by  Kuo  or  Clay  will 

give  good  results  (provided  a  is  correctly  calculated).  However, 

s 

for  backscatter,  or  low  grazing  single  forward  scatter  in  the  case  of 
one  interface,  or  any  type  of  scattering  in  the  case  of  multiple 
layers,  Eq.  (57)  will  not  hold. 

For  the  ocean  bottom  models  to  be  used  in  this  section,  the 
scattered  intensity  is  given  by 


r(  ei,T))R*(>*i,V  )D0*D0 
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where  the  notation  is  identical  to  the  notation  used  in  Eq.  (26),  and 
the  source  strength  has  been  normalized  to  unity.  This  equation  is 
the  extension  to  acoustically  penetrable  rough  surfaces  of  the  potential 
formulation  given  in  the  Final  Report  under  Contract  N00024-69-C-1275. 
The  boundary  value  has  been  modified  by  the  inclusion  of  R(@^,r,), 
the  pressure  density  coefficient  (or  reflection  coefficient).  R(e^, q) 
is  given  by  either  Eq.  (17)  or  Eq.  (19)  according  to  whether  the 
bottom  contains  a  single  or  double  interface . 


Upon  making  the  stationary  phase  calculation  given  in  Section  II, 
Eq.  (58)  becomes 
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The  integrals  given  in  Eq.  (59)  have  been  calculated  in  the  previously 
mentioned  final  report  for  a  surface  with  a  normal  (Gaussian) 
bivariate  height  distribution  and  Gaussian  correlation  function,  and 
they  included  a  realistic  insonification  function  and  the  Fresnel 
phase  approximation. 


The  intensity  scattered  by  an  acoustically  penetrable  rough 
surface  is  then  given  by  [recall  F(£l,#r)  is  defined  by  Eq.  (56)] 
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where 


I2  2V2 

g  =  k  7  h  , 

A  is  the  insonified  area  defined  by  the  -3  dB  down  points. 


M  = 


N  = 
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K  - - go -  is  a  beam  function  constant, 

L  is  the  correlation  length  of  the  upper  surface, 

a  =  cos©,  -  cos©  , 
i  r* 

R,  =  2r  r_.  /(r  +r..  ), 

1  oo  lc/  '  oo  lo/J 

R0  =  2r  r.  /(r  sin2©  +  r,  sin2©.),  and 

2  oo  1&  v  oo  r  lo  l '  ’ 

a  and  £  are  the  semimajor  and  semiminor  axes,  respectively, 
of  the  elliptical  insonified  area,  A  =  jtorp. 


Since  the  important  features  of  Eq.  (6o)  are  rather  difficult  to 
discover  in  its  present  form,  let  us  investigate  these  formulas  for  the 
case  where  the  receiver  is  in  the  specular  direction.  Using  the  defini¬ 
tions  of  M,  N  and  a,  and  some  algebraic  manipulation  gives 
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One  significant  but  often  misunderstood  aspect  of  reflection  and 
scattering  phenomena,  demonstrated  in  Eq.  (6l),  is  the  range  dependence. 
The  n=o  term  is  just  the  coherently  scattered  intensity;  since  this 
term  can  be  considered  to  be  the  reflection  by  a  lossy  plane,  the 
geometrical  acoustics  (optics)  range  dependence,  (r  +r^Q)  >  is 
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obtained  as  expected.  For  n=l,  .  .  .  ,  °°,  the  terms  in  Eq.  (6l) 
represent  the  incoherently  scattered  intensity.  Clearly,  for  N  and 
M  very  large  (recall  that  M  and  N  depend  on  R^  and  R,  respectively), 
rhe  range  dependence  of  the  incoherently  scattered  intensity  will 
also  be  given  by  geometrical  acoustics.  The  condition  of  M  and  N 
large  is  obtained  when  the  insonified  area  is  many  wavelengths  in 
dimension  and  the  roughness  is  not  too  great.  When  M  and  N  are 
small,  Eq.  (60)  becomes 
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Here,  the  incoherently  scattered  pressure  is  now  area  dependent  and 
has  Fraunhofer  range  dependence,  (r00rl0)"  •  The  condition  of  M  and 
N  small  which  leads  to  the  above  equation  is  obtained  when  the  insoni¬ 
fied  area  is  small  with  respect  to  a  wavelength  and/or  when  the 
roughness  features  are  quite  large.  Even  for  a  large  insonified  area 
it  is  interesting  to  note  that  Eq.  (62)  will  be  obtained  for  quite 
rough  surfaces.  The  physical  interpretation  of  this  result  is  that 
for  a  fairly  rough  surface  the  individual  surface  features  behave  as 
independent  soatterers  with  dimensions  on  the  order  of  a  wavelength. 
Equation  (62)  is  also  obtained  if  the  insonified  area,  is  only  a  few 
wavelengths  in  diameter. 

Equation  (6o)  thus  has  two  important  consequences.  In  addition 

to  its  ability  to  predict  the  correct  shape  of  the  forward  scattered 

field,  it  also  predicts  that  the  range  behavior  of  the  scattered 

field  makes  a  transition  from  (r  +r,  )  for  a  very  smooth  surface 
2  '  oo  lo'  J 

to  (r  r  )  Tor  a  very  rough  surface  or  for  a  very  small  insonified 
area.  The  transition  in  the  range  dependence  is  obtained  strictly  as 
a  result  of  using  the  Fresnel  phase  approximation.  This  situation 
should  be  contrasted  with  the  range  dependence  by  all  of  the  other 
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scattering  theories  (see,  for  example.  Refs.  5,6,7,8,9,10,11,12,13,1^, 

16,17)  which  use  only  the  Fraunhofer  phase  approximation,  and  conse- 

-2 

quently  obtain  only  a  range  dependence  of  (r00r]_0)~  •  Obviously,  in 
the  limit  as  the  surface  approaches  a  plane  this  result  will  fail 
(unless,  it  can  be  shown  that  a  real  source  at  infinity  still  insoni- 
fies  only  an  area  of  a  few  wavelengths). 

To  obtain  the  scattering  coefficient  generally  used  in  field 
work,  the  source  strength  and  the  total  travel  path  are  usr  /'  . r 
divided  out.  Since  the  source  strength  has  already  been  nor^lized 
to  unity,  the  scattering  coefficient  of  an  acoustically  penetrable, 
rough  surface  is  given  by 

/en  =  /i  \  (r  +r  )^  ,  (63) 

S  \  8/  OO  lo'  ’  '  ' 

where  <^1^  is  given  by  Eq.  (60).  This  result  will  be  used  in  the 
next  subsection  in  conjunction  with  the  parameter  values  listed  in 
Table  I  to  demonstrate  the  effects  of  an  acoustically  penetrable 
bottom  on  bottom  bounce  sonar  operation. 

C.  Discussion  of  Results 

The  acoustic  parameters  given  in  Table  I  will  be  used  in 
Eq.  (63)  to  study  the  effects  of  some  typical  rough,  acoustically 
penetrable  bottoms  on  bottom  bounce  sonar  operation.  The  only 
additional  quantities  which  must  be  supplied  are  the  rms  heights  and 
slopes  which  might  be  encountered  in  the  different  physiographic 
provinces .  In  the  interests  of  brevity,  only  two  bottom  provinces 
are  considered.  On  a  shallow  bottom  (continental  margin),  greater 
roughness  can  be  expected  than  on  a  deep  bottom  (ocean  basin)  due 
to  the  presence  of  turbidity  currents  and  much  greater  depositional 
activity.  Due  to  the  great  variability  to  be  found  in  the  bottom 
structure,  it  will  be  possible  to  investigate  the  scattering  for 
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only  a  few  typical  but  arbitrarily  chosen  roughness  values.  The 
values  of  the  rms  height  and  slope  which  wi^l  be  used  here  are  given 
in  Table  II.  For  the  single  interface  case,  the  values  of  a  and  s 
represent  the  rough  pax’ameters  at  the  water- sediment  interface.  For 
the  double  interface  case  the  subscript  1  indicates  the  water- sediment 
roughness  parameters,  while  tne  subscript  2  indicates  the  sub-bottom 
roughness  parameters .  Notice  that  the  sub-bottom  roughness  is  here 
the  same  for  both  deep  and  shallow  bottoms . 

In  Dwg.  AS- 71- 389  the  effect  of  roughness  alone  on  a  sonar  beam 
reflected  from  an  acoustically  impenetrable  bottom  in  the  plane  of 
incidence  is  illustrated.  The  sonar  has  a  conical  beamwidth  of 
8  deg  measured  with  respect  to  the  half- power  points  and  is  operating 
at  a  frequency  of  3.5  kHz  (X  a  18  in.). 

TABLE  II 

SOME  TYPICAL  ROUGHNESS  PARAMETERS  FOR 
REALISTIC  OCEAN  BOTTOMS 


Bottom  Province/ 

Model 

RMS  Height, 

a 

RMS  Slope, 
s  =  a/L 

1.  Continental  Margin, 
Single  Interface 

ct  =  12  in. 

s  =  0.1 

2.  Continental  Margin, 
Double  Interface 

=  12  in. 

^2  =  9  in. 

S1  =  0.1 
s2  =  0.05 

3.  Ocean  Basin, 

Single  Interface 

a  =  1  in. 

s  =  0.0085 

*■  *  .  '  ’-i" — — 

4.  Ocean  Basin, 

Double  Interface 

a^=l  in. 

CT2  =  9  in. 

sx  =  0.0085 
s2  =  0.05 
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The  beam  reflected  by  a  perfectly  reflecting  plane  bottom  is 
included  for  comparison.  The  beamwidth  of  the  beam  reflected  by  a 
plane  appears  to  be  broader  than  8  deg  for  the  simple  reason  that 
the  angles  shown  in  the  drawing  are  measured  with  respect  to  the 
point  where  the  beam  axis  intersects  the  surface.  This  makes  the 
beam  appear  to  be  about  twice  as  broad  as  it  would  be  if  the  angles 
were  measured  with  respect  to  the  image  point  where  the  source  can  be 
thought  of  as  actually  being  located.  The  slight  roughness  assumed 
for  a  deep  bottom  has  resulted  in  a  beamwidth  increase  from  8  deg 
to  9  deg  at  30  deg  grazing  incidence  with  no  other  perceptible  effects. 
The  greater  roughness  assumed  for  the  shallow  bottom  has  resulted  in 
a  beamwidth  increase  from  8  deg  to  12  deg  at  30  deg  grazing  incidence 
and  the  loss  of  about  4  d£  from  the  specular  portion  of  the  beam.  It 
should  be  emphasized  that  the  beam  broadening  is  given  with  respect 
to  the  source,  i.e.,  the  ship;  it  appears  to  be  larger  in  Dwg.  AS-71-389 
because  of  the  way  the  angles  are  referenced. 


The  curves  given  in  Dwg.  AS-71-390  represent  the  effect  of  the 
single  and  double  interface  shallow  bottom  model  on  a  bottom  reflected 
sonar  beam  at  15  deg  grazing  angle .  The  beam  which  would  be  reflected 
by  an  acoustically  impenetrable  plane  is  also  included  for  reference. 
For  the  single  interface  case  the  peak  value  of  the  beam  is  displaced 
downward  by  about  3  deg  from  the  beam  axis.  The  reason  for  this 
angular  displacement  of  the  beam  maximum  is  that  the  portions  of  the 
beam  which  are  incident  at  higher  grazing  angles  have  greater  pene¬ 
tration  into  the  bottom  if  >  c,,.,™™.  For  a  low  speed  bottom 

the  peak  would  be  displaced  to  the  opposite  side  of  the  beam  axis  or 
the  beam  would  have  a  null  at  the  angle  of  intromission.  There  is 
also  a  4  dB  loss  in  the  specular  direction  due  to  penetration,  in 
addition  tc  the  k  dB  scattering  loss.  When  an  intermediate  deposi- 
tional  layer  is  present  that  provides  a  good  impedance  match  between 
the  water  and  bottom,  then  the  penetration  losses  increase  drastically. 
Ordinarily,  for  a  smooth  sub-bottom  and  little  absorption  the  peaks 
of  the  multiple  lobes  due  to  the  sub-bottom  interference  would  be 
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very  close  to  the  level  of  the  single  Interface  scattering.  However, 
since  the  sub-bottom  can  also  scatter  acoustic  waves,  the  sub-bottom 
is  a  source  of  scattering  interference  with  the  waves  reflected  by 
the  bottom  itself.  The  amount  of  lobing  that  occurs  depends  on  how 
many  wavelengths  the  sub-bottom  lies  below  the  bottom.  For  the 
assumed  shallow  bottom  model,  the  sub-bottom  lies  about  7  wavelengths 
below  the  bottom,  and  consequently  there  is  considerable  lobing  due 
to  interference.  Drawing  AS-71-391  shows  the  effect  of  the  shallow 
(rough)  bottom  on  a  sonar  beam  incident  at  30  deg  grazing. 

Drawings  AS- 71- 392  and  AS-71-393  show  the  effect  of  the  typical 
deep  ocean  bottom  model  on  a  sonar  beam  incident  at  15  deg  and  30  deg 
grazing,  respectively.  Since  the  water- sediment  interface  was 
assumed  to  be  quite  smooth,  little  beam  broadening  has  taken  place. 
Again  the  main  effects  of  the  bottom  at  15  deg  grazing  have  been  the 
displacement  of  the  beam  peak  and  some  penetration  losses.  In  the 
case  of  an  impedance  matching  intermediate  layer,  there  is  some 
additional  penetration  loss  with  very  little  lobing.  The  lobing 
effect  is  much  less  here  since  the  sub-bottom  lies  only  one  and 
one- half  wavelength  below  the  bottom.  At  30  deg  grazing,  the  main 
effect  of  a  smooth,  acoustically  penetrable  bottom  has  been  to  produce 
a  slightly  narrower  beam  with  about  10  dB  loss  near  specular  angles 
and  only  minor  lobing. 

It  must  be  emphasized  that  the  descriptions  of  shallow  and  deep 
bottoms  given  here  are  strictly  arbitrary;  there  is  no  rule  which 
states  that  a  deep  bottom  must  be  smoother  and  have  thinner  sediment 
deposits  than  a  shallow  bottom.  On  the  average  these  descriptions 
might  hold,  and  at  the  very  least  they  can  be  considered  reasonably 
representative  of  ocean  bottoms. 

In  summary,  several  general  comments  can  be  made  on  the  effects 
of  rough,  acoustically  penetrable  bottoms  on  bottom  bounce  sonar 
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PLANE,  IMPENETRABLE  BOTTOM 
DEEP  BOTTOM,  SINGLE  INTERFA< 


PLANE,  IMPENETRABLE  BOTTOM 
DEEP  BOTTOM,  SINGLE  INTERS CE 
DEEP  BOTTOM,  DOUBLE  INTERFACE 


EFFECT  OF  A  ROUGH,  ACOUSTICALLY  PENETRABLE  BOTTOM  ON 
A  BOTTOM  BOUNCE  SONAR  BEAM  AT  30  deg  GRAZING  ANGLE 


operation.  If  the  bottom  is  rough,  then  directional  ambiguity  is 
introduced  into  the  target  location  due  to  the  beam  broadening 
caused  by  scattering.  Further,  if  the  bottom  is  both  rough  and 
penetrable,  then  considerable  beam  distortion,  or  lobing,  and  level 
decrease  due  to  penetration  may  occur  on  bouncing  a  signal  from  the 
bottom.  Also,  same  of  the  secondary  effects  would  become  more 
important.  For  example,  volume  reverberation  would  increase  in 
duration  and  possibly  level  since  scattering  would  increase  the 
insonified  ocean  volume.  In  addition,  broad  scattered  beams  would 
contribute  to  mutual  interference  among  several  sonar  systems 
operating  in  the  same  vicinity.  Finally,  a  rough  bottom  introduces 
fluctuations  Into  the  scattered  sound  field  as  the  insonified  area 
moves  along  the  bottom.  Of  course,  the  fluctuations  depend  on  the 
density  of  scatterers  in  the  insonified  area.  If  there  are  very 
many  scatterers  in  the  insonified  area,  then  the  fluctuations  will 
be  small  as  the  insonified  area  moves  along  the  bottom,  since  each 
return  will  represent  a  large  statistical  sample.  For  this  reason 
the  fluctuations  in  a  Bignal  scattered  by  a  deep  bottom  will  be  less 
than  in  a  signal  scattered  by  the  identical  bottom  at  a  shallower 
depth,  since  the  insonified  area  increases  with  depth. 
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IV.  THE  IN  SITU  DETERMINATION  OF  THE  SUITABILITY  OF  THE 
OCEAN  BOTTOM  FOR  BOTTOM  BOUNCE  SONAJR  OPERATION 

There  are  many  operational  circumstances  in  which  one  would  like 
to  use  the  bottom  bounce  mode,  but  knowledge  of  the  suitability  of 
the  bottom  for  such  sonar  operation  is  lacking.  Regional  acoustic 
atlases  have  been  compiled,  but  at  best  these  can  supply  only  very 
broad  outlines  as  to  the  usefulness  of  the  bottom  bounce  mode  at  a 
given  location.  Ideally,  what  a  shipboard  sonar  operator  desires 
would  be  a  very  simple  method  to  determine  the  characteristics  of 
the  bottom  using  readily  available  equipment. 

If  there  is  only  one  ship  available,  then  the  bottom  measurements 
must  be  made  in  the  backscattering  configuration.  The  question  then 
naturally  reduces  to:  What  can  be  predicted  about  the  forward  scatter 
behavior  from  measurements  of  the  backscattering?  To  be  able  to  answer 
this  question  requires  that  the  theoretical  foundations  of  forward  and 
backscattering  be  well  understood. 

The  theoretical  formulation  of  scattering  developed  at 

Applied  Research  Laboratories  (ARL)  predicts  forward  scattering 

accurately  if  the  rms  height  and  slope  of  the  scattering  surface  are 

known.  Also,  backscattering  predictions  are  quite  good  from  normal 

incidence  down  to  50  deg  grazing  angles.  Below  50  deg,  the  typical 

flattening  out  of  the  backscattering  curves,  which  is  observed  in 

all  backscatter  measurements  (acoustic,  radar,  and  laser),  is  not 

predicted  by  the  usual  theories.  It  was  first  postulated  by 
20 

Kur'y&nov  7  that  the  correct  agreement  between  measured  and  predicted 
backscattering  could  be  obtained  if  it  were  assumed  that  the  surface 


was  composite  in  nature;  that  is,  in  addition  to  the  roughness  features 

which  are  easily  measured,  there  occurs  on  most  natural  surfaces  a 

small  scale  roughness  which  is  usually  characterized  by  a  much  larger 

rms  slope  than  that  of  the  large  scale  roughness.  The  concept  of 

scattering  due  to  large  scale  and  small  scale  roughness  features  has 

received  additional  theoretical  and  experimental  support  in  the  papers 

of  Beckmann,^0  Fuks,"^  Fung  and  Chan,^  Volovova  and  Zhitkovskiy, 
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Zhitkovskiy  and  Lysanov,  and  Schmidt.  Some  natural  examples  for 
this  type  of  feature  would  be  ripple  on  swell,  boulders  strewn  on  a 
hill,  etc.  The  first  theories  developed  using  the  composite  surface 
concept  gave  good  predictions  of  backscattering  after  all  of  the 
parameters  had  been  adjusted,  but  in  many  cases  for  the  same  parameter 
values  only  nonsense  could  be  obtained  for  forward  scattering.  The 
reason  for  these  failures  can  be  traced  directly  to  the  approximations 
made  in  the  treatment  of  the  surface  slopes. 

The  composite  theory  developed  at  ABL  utilizes  an  exact  slope 
treatment  and  can,  consequently,  predict  both  forward  and  backscatter¬ 
ing  from  the  same  set  of  parameters.  Unfortunately,  the  small  scale 
features  usually  cannot  be  measured;  hence,  they  can  only  be  determined 
by  adjusting  the  parameters  until  the  theory  fits  experiment.  This 
requires  a  certain  amount  of  faith  on  the  part  of  the  theoretician. 

It  should  be  noted  that  this  procedure  is  not  as  weak  as  it  sounds. 
Since  the  large  scale  roughness  is  well  known,  only  the  small  scale 
parameters  must  be  guessed.  However,  the  guess  obtained  for  one 
surface  at  one  frequency  gives  a  good  fit  for  all  the  other  surfaces 
and  frequencies,  which  instills  considerable  faith  in  the  procedure. 

The  most  surprising  aspect  of  the  small  scale  roughness  features 
is  that  their  effect  is  only  evident  in  low  grazing  angle  back- 
scattering;  forward  scattering  and  high  angle  backscatter  are  largely 
determined  by  the  large  scale  roughness  features. 
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The  backscattering  by  an  impenetrable  surface,  which  has  small 
scale  roughness  (x/a  »  100 )  with  large  rras  slopes,  is  shown  in 
Dwg.  AS- 71- 555.  Notice  the  almost  constant  level  from  75  deg  on 
down  to  20  deg  grazing.  The  backscattering  level  seen  in  this  curve 
can  be  adjusted  to  any  arbitrary  level  by  merely  changing  the  values 
of  the  rms  height  and/or  slope .  Examination  of  the  forward  scatter¬ 
ing  curve  for  the  same  parameter  values  reveals  that  the  relatively 
large  change  in  the  distribution  of  energy  seen  in  the  backscatter 
direction  is  obtained  at  a  cost  of  only  l/3  dB  in  the  forward  specular 
direction.  While  the  curve  given  in  Dwg.  AS-71-355  was  derived  from 

theoretical  considerations,  experimental  curves  of  like  shape  have 

35 

been  obtained  by  Nolle  et  al.  for  a  very  similar  experimental 
situation.  No  attempt  was  made  in  Dwg.  AS-71-355  to  match  the  physical 
parameters  of  Nolle' s  experiments. 

If  the  backscattering  for  impenetrable  surfaces  composed  of  the 
single  small  scale  roughness  but  various  large  scale  roughnesses  is 
plotted,  one  obtains  the  curves  shown  in  Dwg.  AS-71-556.  If  the  small 
scale  parameters  are  adjusted  slightly,  but  the  large  scale  parameters 
held  fixed,  the  curves  shown  in  Dwg.  AS- 71-357  are  obtained.  The  only 
change  seen  is  that  the  limiting  level  at  low  gracing  angles  has  been 
adjusted  upward  5  dB. 

The  significance  of  Dwg.  AS-71-357  is  that  at  grazing  angles 
less  than  45  deg  no  information  on  the  large  scale  roughness  features 
is  retained,  since  the  backscattering  at  those  angles  is  controlled  by 
small  scale  roughness.  Consequently,  if  only  the  backscattering  data 
for  a  single  frequency  and  at  grazing  angles  less  than  45  deg  were 
available,  no  concisions  about  the  large  scale  roughness  could  be 
reached,  and  hence  no  prediction  of  forward  scatter  could  be  made. 
(Remember,  it  is  only  the  large  scale  roughness  which  affects 
forward  or  specular  scattering.)  In  Dwg.  AS-71-356,  it  is  seen  that 
no  information  on  the  large  scale  roughness  is  retained  at  grazing 
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angles  less  than  30  deg.  it  should  be  noted  that  the  scattering  due 
to  the  small  scale  roughness  is  frequency  dependent,  in  some  cases 
varying  as  frequency  to  the  fourth  power;  therefore,  the  masking 
effect  of  the  low  grazing  angle  backscatter  may  not  be  as  pronounced 
at  other  (principally  lower)  frequencies.  The  effects  due  to  the  large 
scale  features  appear  to  be  largely  frequency  independent*  ’  hence, 
the  information  about  the  large  scale  features  will  be  retained  if 
frequency  adjustments  are  made  to  reduce  the  small  scale  scattering. 

In  view  of  the  difference  in  the  levels  at  normal  incidence  for 
the  different  roughnesses  in  Dwg.  AS- 71- 356  and  AS- 71- 357-  the  obvious 
question  might  then  be  asked:  Suppose  the  data  for  normal  incidence 
(90  deg  grazing)  were  available,  could  anything  be  inferred  from  these 
data?  First,  recall  that  the  surfaces  considered  in  these  figures  are 
acoustically  impenetrable.  For  an  ocean  bottom  the  return  at  normal 
incidence  will  depend  on  both  the  penetrability  and  the  roughness  of 
the  surface  as  given  in  Eqs.  (6l)  and  (63).  If  the  acoustic  parameters 
of  the  bottom  are  unknown,  then  it  is  impossible  to  determine  the 
relative  influence  of  the  topography  and  the  acoustic  penetrability  on 
the  scattering.  Consequently,  it  will  not  be  known  if  the  loss 
measured  at  normal  incidence  was  due  primarily  to  scattering  or  to 
penetration.  The  importance  of  this  point  was  discussed  in  the 
example  given  in  the  introduction  to  Chapter  III. 

It  was  pointed  out  in  Chapter  II,  in  the  discussion  of  the 
stationary  phase  evaluation  of  R (0^),  that  for  backscatteriog, 

R(0i)  is  given  by  a  constant  ( |R(0  )|i  l)  for  the  single  interface 
case.  Consequently,  the  shape  of  the  backscattering  curve  from  a 
penetrable  surface  remains  constant,  with  only  the  level  shifted 
downward.  Therefore,  for  a  fixed  frequency,  it  is  necessary  to  have 
the  backscattering  curve  for  the  entire  range  of  grazing  angles 
from  90  deg  to  10  deg  (or,  at  least  from  90  deg  to  40  deg)  if 
any  inferences  about  the  large  scale  roughness  are  to  be 
drawn.  The  implication  of  these  facts  is  rather  severe  if 
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it  was  hoped  that  forward  scattering  performance  of  a  bottom  bounce 
sonar  could  be  predicted  from  simple  observations  of  the  backscatter- 
ing  return  level  from  the  sonar  and  the  on-board  depth  sounder  since 
these  pieces  of  equipment  will  not  cover  the  requisite  angles  or 
frequencies . 

At  present,  it  seems  that  the  only  way  to  obtain  data  on  the 
bottom  roughness  from  backscatter  measurements  using  present  ships' 
equipment  will  require  some  fairly  sophisticated  signal  processing. 
Even  this  method  will  require  that  the  return  from  a  very  small 
region  of  the  bottom  (fractions  of  a  correlation  length)  must  be 
resolved  by  the  receiver. 

If  a  wideband  source,  for  example  an  explosive  charge,  is  used, 
then  the  ability  to  differentiate  between  the  bottom  losses  arising 
because  of  scattering  and  acoustic  penetration  increases.  This 
technique  is  complicated  by  the  necessity  of  using  wideband, 
multichannel  filtering  and  recording  equipment. 

The  separation  of  the  losses  due  to  scattering  and  bottom 
penetration  is  based  on  a  presumed  difference  in  the  frequency 
dependence  of  the  two  effects.  The  scattering  is  assumed  to  be 
frequency  dependent  while  the  acoustic  penetration  is  not.  It  is 
well  known  experimentally  and  theoretically  that  if  the  bottom 
features  are  large  with  respect  to  the  wavelength,  then  the  scatter¬ 
ing  is  frequency  independent  and  varies  only  as  the  slopes  of  the 
scatterer  and  not  the  heights.  This  effect  is  associated  with  the 
high  frequency  (in  reality,  large  g(g  -  k  7  h  )]  limit  of  scattering. 
To  obtain  roughness  information  from  a  wideband  source  requires  that 
the  frequencies  present  should  represent  both  the  high  frequency 
limit  and  the  low  frequency  limit.  The  data  recorded  using  explo¬ 
sive  sources  are  usually  contained  in  the  5  octave  band  from  1  to 
30  kHz.  These  frequencies  represent  wavelengths  of  2  to  60  in. 
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Based  on  the  rather  Incomplete  knowledge  of  the  roughnesses  to  be  found 

in  the  ocean  (see  Chapter  III),  the  frequency  band  of  1  to  30  kHz 

adequately  represents  the  high  frequency  limit  but  is  only  marginally 

close  to  a  low  frequency  limit.  However,  for  many  ocean  regions  the 

1  to  30  kHz  band  provides  sufficient  frequency  dependence  to  determine 

32  3k 

the  relative  effects  of  roughness.  ’ 

The  assumption  that  the  acoustic  penetration  of  the  bottom  is 
frequency  independent  holds  only  for  a  bottom  from  which  the  return 
comes  primarily  from  a  single  interface.  If  the  bottom  consists  of 
several  layers  with  different  acoustical  properties,  the  angular 
dependence  of  the  reflection  coefficient  will  be  complex  and  in  every 
instance  will  be  associated  with  the  parameters  of  the  layers  and  the 
frequency  of  the  sound  wave.  Similarly,  if  the  absorption  is  strongly 
frequency  dependent,  interpretation  of  the  bottom  returns  is  again 
rendered  difficult . 

Considerable  qualitative  data  are,  however,  provided  by  the 

short  pulses  which  are  associated  with  explosive  sources.  Zhitkovskiy"^ 

describes  a  method  of  measuring  the  distances  between  irregularities 

37 

on  the  ocean  floor.  Brekhovskikh  shows  data  from  which  much 
qualitative  infonnation  on  the  character  of  the  bottom  can  be  deter¬ 
mined  by  a  fairly  inexperienced  operator  after  only  a  cursory 
inspection. 

In  summary,  the  feasibility  of  _in  situ  determination  of  ocean 
bottom  characteristics  depends  on  the  available  measurement  apparatus. 
The  bottom  characteristics  of  interest  have  been  shown  to  be  the  large 
scale  topography  of  the  bottom  and  the  acoustic  penetrability  (which 
depends  on  the  physical  composition  and  structure  of  the  bottom) .  To 
adequately  determine  these  characteristics  from  simple  measurements 
of  the  backscattered  signal  level  requires  at  least  the  range  of 


83 


grazing  angles  from  90  deg  to  40  deg.  The  confidence  in  this 
determination  will  be  substantially  increased  if  the  data  for  several 
different  frequencies  are  also  available.  The  equipment  presently 
available  in  tne  form  of  on-board  sonars  and  depth  sounders  is  not 
adequate  for  this  determination.  If  additional  equipment  in  the  form 
of  wideband,  multichannel  filtering  and  recording  equipment  is 
available,  then  explosive  sources  may  be  used  to  obtain  the  necessary 
angular  and  frequency  measurements.  However,  in  many  cases,  the 
reduction  of  these  raw  data  and  the  subsequent  prediction  of  forward 
scattering  require  sophisticated  analysis  by  experienced  personnel. 

The  ideal  situation  would  be  to  have  a  profilometer  which  has  a 
sufficiently  narrow  beam  so  that  the  insonified  area  of  the  bottom 
would  be  a  few  correlation  lengths .  The  time  record  of  this  return 
would  then  represent  a  bottom  profile  with  sufficient  resolution 
to  determine  the  topographical  features  of  interest  in  forward 
scattering.  Further,  if  the  source  were  calibrated,  then  the  level 
of  the  returned  signal  could  be  used  to  determine  the  acoustic 
penetrability  of  the  bottom.  Under  these  circumstances  the  forward 
scattered  signal  could  be  predicted. 
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V.  ACOUSTIC  FIELD  IN  A  SURFACE  DUCT  WITH  A  ROUGH  SURFACE 

A.  Introduction 

The  acoustic  field  is  derived  for  a  cw  point  source  located  in 

an  inhomogeneous  medium  with  a  rough  surface.  Application  is  made 

to  the  case  of  a  surface  duct,  and  examples  ere  given  for  two  cases 

(bilinear  and  Epstein).  The  effects  of  surface  roughness  are  intro- 

38-  40 

duced  by  means  of  an  effective  reflection  coefficient.  Lapin 

41-43 

and  Lysanov  show  that  a  statistically  rough  boundary  can  be 

approximated  by  a  flat  boundary  with  some  effective  reflection 

coefficient  in  the  case  of  waveguide  propagation,  but  the  form  of 

this  coefficient  will,  in  general,  be  different  from  the  expression 

43 

for  reflection  frcm  a  single  uneven  surface,  lysanov  explicitly 
carries  out  the  calculations  for  the  case  of  a  Gaussian  surface 
overlying  a  bilinear  surface  duct. 

44 

Recently  Bucker  has  given  an  excellent  treatment  of 
propagation  in  a  layered  waveguide  with  lossy  boundaries.  Bucker, 
however,  gives  the  effect  of  roughness  in  terms  of  a  plane  wave 
reflection  coefficient  S,  which  does  not  allow  for  the  variable 
nature  of  the  medium.  The  reflection  is  assumed  to  take  place  from 
a  pseudolsovelocity  layer  which  is  allowed  to  shrink  to  zero.  Bucker 
makes  no  attempt  to  drive  S  from  the  statistical  nature  of  the  surface. 
He  assumes  various  values  for  the  dB  loss  per  surface  reflection  SL 

n 

and  then  calculates  the  reflection  coefficient  S  from  the  equation 

SLn  =  -20  log10  | S j  (64) 
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The  propagation  model  that  Bucker  foimulates  assumes  an  impedance 
boundary  condition.  The  value  for  the  impedance  is  determined  from 
the  plane  wave  reflection  coefficient.  In  order  to  find  initial 
approximate  eigenvalues,  Bucker  relates  each  mode  to  a  ray  and 
arrives  at  a  boundary  attenuation  in  the  form  of  an  imaginary  component 
to  the  eigenvalues.  Specifically,  he  obtains 


SL  +  BL 
n 


n 


H 


n 


=  -20  log10  e 


-k" 
n 


(65) 


where  BL  is  the  bottom  loss  of  the  ray  associated  with  mode  n, 
n 

is  the  ray  horizontal  cycle  distance,  and  k”  is  the  approximate 

imaginary  part  of  the  eigenvalue  which  is  due  to  the  lossy  boundaries. 

4.5 

Bucker  assumes  various  values  for  SL  since  the  values  for  SL  given 
46  n  n 

by  Marsh  seem  to  be  too  high  for  the  experimental  cases  he  treats. 

The  rough  surface  propagation  problem  has  also  been  treated  by 

47 

Clay.  Clay  shows  that  surface  scattering  produces  an  increase  in 

mode  attenuation  and  a  decrease  in  coherence  of  the  signal.  Ray 

48-55 

theory  has  been  applied  to  this  problem  by  various  authors.  A 

useful  survey  of  both  ray  theory  and  wave  approaches  has  been  given 
54 

by  Schulkin.  The  present  report  will  follow  the  treatment  given  by 
Bucker,  but  will  use  Green's  functions  and  attempt  to  arrive  at 
expressions  for  the  surface  loss  per  bounce  which  relate  directly  the 
ray  angle,  frequency,  and  peak  to  trough  heights. 

B.  Green's  Function  Solution 

The  wave  equation  for  a  point  source,  with  angular  frequency  u>, 
located  in  a  layered  inhomogeneous  medium  is  given  by 

V2P  -  =  -io)Q  5(r-r  )  e~mt  ,  (66) 

c2(z)  dt2  °  ° 


06 


where  P  represents  the  pressure,  c(z)  is  the  sound  velocity 
(variable  in  the  coordinate  z),  and  5(r-r  )  is  the  three-dimensional 
Dirac  delta  function.  The  strength  Qq  of  the  source  is  equal  to  the 
density  times  the  volume  velocity,  that  is,  mass  flow  per  second.  If 
cylindrical  coordinates  (r, 0,z)  are  assumed  with  azimuthal  symmetry 
and  a  time  factor  of  exp  (-icut)  is  suppressed,  the  following  equation 
is  obtained: 


1  5 

r  ^r 
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♦  A  ♦  ^)P  . 

8z 
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(67) 


where  P  now  represents  the  time  independent  pressure.  The  wave 
number  k(z)  is  defined  as  a/c(z)  and  the  source  is  located  at  z  =  zq 
and  r  =  o,  as  indicated  in  Dwg.  AS-70-767.  The  vertical  depth 
coordinate  z  varies  from  0  S  z  S  ®,  and  the  range  coordinate  r  varies 
from  0  <  r  <  «.  The  boundary  conditions  imposed  on  Eq.  (67)  are  that 

a)  P  must  satisfy  a  radiation  condition  for  r  -» <*>  and  z  -» °° 


») 


=  7 


z=o 


(68) 


c)  P  and  be  continuous  across  any  discontinuities  in  the 
velocity- depth  profile  (such  as  the  bilinear  profile). 

The  statistical  nature  of  the  surface  at  z  =  o  will  be  introduced  by 
means  of  an  effective  reflection  coefficient  which  then  determines  7. 


In  the  Green's  function  approach,  Eq.  (67)  is  separated  into 
the  following  forms: 
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[Mrfe)  +  Xr]0i<i'’x)=^)  <69) 

+  k2(z)  -  xlo2(z,z0,-k)  =  -±a,q  l(z-z:_j  ,  (70) 


where  X  is  the  separation  constant.  ”>y  the  resolvent  Green's  function 
technique,  the  solution  of  Eq.  (67)  is  given  by  a  complex  convolution 
of  G^  and  G^: 


P(r,: 


'zo} 


2*i  fr  Gl^r> 


X)G0( 


z,z 


•X)dX 


(Tl) 


where  the  contour  C  separates  the  singularities  of  the  Green's  functions. 

The  solutions  of  Eq.  (69)  and  Eq.  (70)  that  satisfy  the  required 
boundary  conditions  are  given  by 


G^r,!)  (  {r)(  $  =  XX/2,o  <  arg  X  <  at) 


(72) 


n 2(z,|)[n  (z  0  -  X"n  (z  ,5)] 

Ga(.,V-X).-l - Srn-nTT-  Zo<Z<“  >  ™ 


.  Mz  ,|)[n  (z,|)  -  X'n  (z,|)] 

GcAz>z~>-X)  =  - T7TZ — r— * -  o  <  z  <  z_ 


where  H^( |r)  is  the  Hankel  function  of  the  first  kind,  and 


nj_(0,  |)  -  7^(0,!) 

n£('°,5)  -  Xn2^0,  |) 


(74) 
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The  Wronskian  of  ng  and  n  is  represented  by  W(n?,n1),  where  n?  is 
a  solution  to  the  homogeneous  part  of  Eq.  (70)  and  has  outgoing  waves 
at  z  and  where  n^  has  outgoing  waves  at  z  -♦-*>.  The  primes 

indicate  differentiation  with  respect  to  z. 

When  Eq.  (72)  and  Eq.  (73)  are  used  in  Eq.  (7l),  the  solution 
for  Eq.  (67)  is  obtained  as 


P(r,z,zQ) 


o£o  f  n2(z<,  5)[ni(z>,|)  -  X‘n2(z>,  |r)^d| 

Tit  J  c  W(n2,ni) 


(75) 


where  z  and  z>  denote  the  smaller  or  larger,  respectively,  of  the 

variables  z  and  z  . 

o 

When  Eq.  (75)  is  integrated  by  Cauchy’s  residue  theorem,  the 
normal  modes  (plus  any  branch  line  integrals)  are  obtained.  The 
poles  yielding  the  normal  modes  are  given  by  the  zeroes  of  the 
expression 


[03(0,5)  -  7n2(0,|)] 


(76) 


The  normal  mode  expansion  can  be  written  as 
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or 


Equation  (77)  can  also  be  written  as 

(79 


In  general  the  eigenvalues  will  be  complex.  The  lossy  boundaries 
add  an  extra  imaginary  component  to  the  complex  eigenvalue  (in 

the  surface  duct  case,  boundaries  mean  the  surface  and  any  velocity- depth 

discontinuities.)  Approximate  values  for  k"  can  be  found  in  the 

.  44  55  n  56  N 

following  manner  (Bucker,  Urick,  and  Watt  ): 

SL  -k" 

s  -20  log1Q  e  n  ,  (80) 
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where  SL  ir  the  surface  loss  per  bounce,  and  H  is  the  ray  horizontal 
n  *  ’  n 

cycle  distant. Q .  SL^  is  defined  by 

SLn=  -20  log10  |S|  (82) 

where  S  is  t*-  nlane  wave  reflection  coefficient  from  a  pseudoisovelocity 
surface  layer  that  is  allowed  to  shrink  to  zero. 


SL 

kn  =  '8T5B6  h” 
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Explicitly,  S  is  related  to  the  boundary  condition  at  the  surface 
through  the  formula 

r  S+l 

2  2  2 

where  (  ^a)  =  (to/vg)  -  and  vg  is  the  sound  velocity  at  the 

surface.  Therefore,  if  either  S  or  7  is  a  known  quantity,  it  is 
possible  to  find  the  attenuation  (k|j)  from  Eq.  (76).  In  this  report 
the  plane  wave  reflection  coefficient  will  be  assumed  to  be  known. 

The  mode  attenuation  factor  k^  is  approximately  given  by 

-SOlo^lSl 
n  8.686  H 

n 

From  classical  scattering  theory  the  coherent  reflection 
coefficient  for  a  Gaussian  surface  is  given  by 


(84) 


(83) 


-2(k  sine  0)^ 

S  =  -e  n  ,  (85) 

where  9  is  defined  by  the  condition  cos  9  =  [v  Re(  f  )/to] ,  and 
a  is  the  rms  height  value  for  the  rough  surface.  If  a  peak  to  trough 
height  h  is  known,  then  a  good  approximation  (Longuet- Higgins)  for  a 
is  given  by  the  equation 


_  h 

2J2 


(86) 


Equation  (82)  will  now  be  calculated  for  one  of  the  cases  Bucker 

considered.  He  found  that  a  reasonable  value  for  SL  was  between  zero 

n 

and  one  dB  per  bounce,  yielding  a  value  of  S  between  -1  and  -0.8913- 
For  a  peak;  to  trough  height  of  2  ft  and  an  average  grazing  angle  of 
3  deg,  Eq.  (85)  yields  a  value  of  S  equal  to  -0.92  for  a  frequency  of 
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1.5  kHz.  Therefore,  Eq.  (85)  appears  to  be  a  valid  method  for 
introducing  the  surface  roughness  into  waveguide  problems. 

46  57 

Marsh  and  Marsh,  Schulkin  and  Kneale  have  arrived  at 

another  expression  for  the  surface  loss  per  bounce  using  the  Neumann- 

Pierson  spectrum.  Their  expression  for  SL^  is  given  by 

SLn  =  -10  log  Q(l  -  0.0013  b^h1/10!)  ,  (87) 

where  h  is  the  average  trough  to  crest  height  in  feet,  L  is  the 
isothermal  layer  depth  in  feet,  and  b  =  fh  is  the  frequency-wave 
height  product  in  kHz  x  ft.  Since  Marsh's  formula  for  SLn  does  not 
depend  on  the  grazing  angle  0^,  Eq.  (87)  is  not  believed  to  be  a 
valid  formula  for  introducing  the  losses  due  to  surface  scattering. 

Computer  programs  for  propagation  loss  versus  range  are 
presently  being  completed  for  two  surface  duct  profiles.  The 
profiles  (Epstein  and  bilinear)  are  shown  in  Dwg.  AS-71-188.  (The 
figure  has  been  taken  with  permission  from  the  paper  by  Bucker  and 

cO 

Morris.  )  Various  surface  conditions  are  being  considered  for 
these  two  velocity- depth  profiles,  and  propagation  loss  curves  will 
be  published  as  soon  as  the  computer  programs  are  completed.  The 
Epstein  profile  is  a  five- parameter  velocity  function  which  has  be<-\ 
examined  for  ray  theory  solutions  and  normal  mode  solutions. 

Explicitly,  the  Epstein  profile  is  given  by 


where  the  quantities  A,  B,  D,  z^,  and  K  are  parameters . 
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The  bilinear  profile  is  also  a  well  known,  extensively  treated, 
50 

velocity  function.  It  is  given  by 


c 


In  Dwg.  AS- 71- 188  the  parameters  are  given  by  c  =  k9^0.9^  ft/sec, 

z  =  326.67  ft,  7  =  0.0144,  and  71  =  -0.k38.  These  values  were 

a  0  1  50 

taken  directly  from  Pedersen's  article. 

Table  II  Illustrates  the  effects  of  surface  roughness  on 
the  first  three  modes  of  the  Epstein  surface  duct  considered  by 
Bucker  and  Morris.  The  values  for  Re(c  )  In  column  two  do  not 
correspond  to  the  roots  obtained  from 

Sq)  -  ^  =  0  » 

but  correspond  to  the  pressure  release  condition 

= 0  • 

However,  for  the  first  few  modes  the  values  of  Re(cn)  will  nearly 
be  the  same.  In  column  three  the  peak  to  trough  heights  are  given 
from  zero  to  6  yd.  The  values  for  the  rms  heights,  calculated  from 
Eq.  (86),  then  range  from  zero  to  2.12  yd.  Column  four  gives  the 
surface  loss  per  bounce  as  found  from  Eq.  (82).  The  grazing  angles 
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SURFACE  LOSSES  FOR  FIRST  THREE  MODES 


for  the  equivalent  rays  are  given  in  column  five.  In  column  six  the 
reflection  coefficients  as  calculated  from  Eq.  (85)  are  given. 

Column  seven  gives  the  horizontal  range  values  in  kiloyards  as  calcu¬ 
lated  from  the  following  expression  given  by  Pedersen  and  White: 

Hn  =  Ha"1/2  in  |d|  -  in  ^(a^)1/2  +  ajj 


where  all  quantities  on  the  right  of  the  equation  are  defined  in 

Ref.  oO,  and  where  the  subscripted  quantities  are  evaluated  at 

the  surface.  The  mode  attenuation  factor  k"  is  calculated  from 

n 

Eq.  (84)  and  is  given  in  column  eight.  The  reason  for  the  low  losses 
due  to  surface  roughness  is  because  of  the  small  grazing  angles  and 
the  fairly  low  frequency. 


97 


APPENDIX  A 

MODELING  A  PENETRABLE  SURFACE:  SELECTION 
OF  MATERIALS  AND  PARAMETER  MEASUREMENT 

fey 

Helge  Wieder*  and  Pat  Welton 

The  use  of  model  surfaces  in  experimental  scattering  measurements 
is  a  well  established  procedure.  The  models  heretofore  used  in  scatter¬ 
ing  work  at  Applied  Research  Laboratories  have  been  made  of  pressure 
release  materials.  These  scattering  surfaces  are  representative  of 
the  ocean  surface  or  an  extremely  gaseous  ocean  bottom.  The  scattering 
by  the  pressure  release  model  surfaces  is  influenced  only  by  the 
topography  of  the  surface.  The  roughness  of  the  topography  of  the  models 
varies  from  quite  smooth  to  very  rough,  and  by  using  different  frequencies 
the  models  can  be  scaled  to  represent  regions  of  the  ocean  bottom  ranging 
from  abyssal  plain  areas  to  the  midocean  ridges. 

In  modeling  acoustically  penetrable  ocean  bottoms,  not  only  should 
the  expected  topography  be  well  modeled,  but  also  the  physical  parameters 
and  structure  of  the  model  should  bear  a  reasonably  direct  relationship 
to  realistic  ocean  values .  While  it  is  easy  to  obtain  the  topography 
desired  by  simply  making  molds  of  the  pressure  release  models,  the 
scaling  of  the  physical  parameters  required  considerable  effort  in  the 
selection  and  measurement  of  the  various  materials  available.  In  this 
appendix  the  selection  criteria  of  a  modeling  material  and  the  actual 
measurement  of  the  physical  parameters  of  this  material  will  be 
discussed. 

_ 

Visiting  Exchange  Scientist  from  the  Federal  Republic  of  Germany 
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A. 


Selection  Criteria 


Recall  that  the  experimental  work  which  is  to  be  performed  under 
this  program  is  the  measurement  of  scattering  by  an  acoustically 
penetrable  liquid  bottom  of  known  statistical  description  as  a  function 
of  angle  and  frequency.  If  the  statistical  parameters  are  to  remain 
constant,  then  the  model  material  must  be  permanently  moldable.  This 
requirement  eliminates  the  use  of  the  actual  muds  and  unconsolidated 
sediments  found  on  the  ocean  bottom  in  the  making  of  the  models,  since 
the  statistics  would  be  unstable.  When  considering  a  permanently 
moldable  material  the  first  criteria  which  must  be  applied  is  the 
approximation  of  this  material  to  a  liquid.  Essentially,  this  implies 
that  the  shear  velocity  must  be  small  with  respect  to  the  compressions! 
velocity.  Some  of  the  first  materials  which  come  to  mind  are  trans¬ 
ducer  window  materials.  The  agents  are  polyurethane  compounds  which 
behave  like  rubber  in  that  they  will  only  support  a  very  small  shear 
wave  component,  and  hence  very  little  conversion  from  compressional  to 
shear  mode  occurs.  These  materials  have  the  further  advantage  of  having 
a  pc  product  ( p  is  the  density,  and  c  is  the  speed  of  sound  in  the 
material)  which  is  very  close  to  that  of  water  and  which  allows  these 
materials  to  be  used  as  a  binding  agent  for  sedimentary  materials.  In 
this  manner  a  moldable  compound  can  be  made  which  has  pc  values  veiy 
close  to  those  actually  encountered  in  the  ocean  bottom  but  for  which 
the  binding  agent  is  practically  acoustically  invisible. 

The  binding  agent  is  Scotchcast  221  polyurethane,  while  tne 
filler  material  is  No,  5  sandblasting  sand.  The  acoustic  characteristics 
of  the  penetrable  material  can  be  altered  by  varying  the  proportion  of 
filler  material  to  binding  agent.  Number  5  sandblasting  sand  was 
selected  as  a  filler  material  because  it  is  the  finest  sand  commercially 
available.  A  very  fine  sand  is  necessary  to  prevent  settling  of  the 
filler  during  the  curing  of  the  binding  agent.  The  only  disadvantage 
of  using  sand  as  a  filler  material  is  that  a  pc  product  less  than  that 
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of  water  cannot  be  obtained.  However,  this  defect  can  be  remedied  by 
using  the  very  fine,  commercially  available  plastic  beads  which  are 
used  to  make  reflective  paint. 

It  should  be  noted  that  as  the  proportion  of  filler  material 
to  binding  agent  is  increased,  the  rigidity  and  hence  the  shear 
velocity  of  the  compound  increases,  and  the  approximation  to  a  liquid 
bottom  is  diminished. 

B.  Measurement  of  the  Acoustic  Parameters 


Appropriate  samples  were  made  for  each  of  the  separate  parameter 
measurements.  Each  sample  was  placed  in  a  vacuum  chamber  prior  to 
curing  to  remove  any  air  bubbles  entrained  during  the  thorough  mixing 
of  fillei^aterial  and  binding  agent.  The  measurements  described  in 
this  section  were  made  for  pure  Scotchcast  221  or  for  Scotchcast  and 
sand.  No  samples  of  Scotchcast  and  plastic,  beads  nave  been  prepared 
as  yet. 

1 .  Shear  Velocity  in  Pure  Scot-hcust 

The  shear  velocity  determined  oy  measuring  the  Lame 
constant,  g,  and  by  using  the  equation 


The  value  of  4  was  determined  according  to  two  standard 
mechanical  tests.  'The  first  determination  was  made  by  the  method  of 
the  torsional  pend *31  urn. 

A  cylindrical  sample  of  Scotchcast  221  of  length  8  in.  and 
radius  l/8  -in.  was  rigidly  clamped  at  the  upper  end,  while  a  metal 
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weight  of  moment  of  inertia  6  was  rigidly  fixed  at  the  lower  end.  The 
metal  weight  was  displaced  a  small  angular  distance  from  equilibrium 
and  was  released  so  that  it  described  torsional  oscillations.  The 
period  of  the  oscillations  was  measured.  The  value  cf  p  is  oh  Gained 
from 


^  - 


to 


(A2) 


where  u>  is  the  angular  frequency  of  the  oscillations,  l  and  r  are  the 
length  and  radius  of  the  sample,  and  Q  is  the  moment  of  inertia  of 
the  suspended  weight,  which  can  be  easily  determined  from  standard 
formulas . 

The  second  method  of  determining  p  was  a  standard  deformation 
test.  A  cylindrical,  sample  of  Scotchcast  (diameter  3  in.,  length 
4.1?  in.)  was  deformed  in  a  hydraulic  press  by  a  force  in  the  direction 
of  the  axis.  Length  and  diameter  of  the  sample  were  measured  as  a 
function  of  force.  From  these  measurements  Young's  modulus  E  and 
Poisson's  ratio  a  were  determined.  Then 

=  2(l+a)  * 

The  values  of  4  determined  from  the  two  methods  were  very  close: 

6  N 

PT  =  1.53  x  10  ~ 

m 

Hjj  -  1.31  Xio6^ 
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then  the  shear  velocity 


A 

if  the  density  of  Scotchcast  221  is  1.06  g/cm , 
is 

cg  =  35.3  m/sec  =*  115  ft/sec 

Subsequent  error  analysis  indicated  the  values  of  u  were  probably 
accurate  to  within  ±2$. 

2.  Density  and  Porosity  Measurements  of  Mixtures  of 

Scotchcast  221  and  Sand 

The  No.  5  sand  used  in  these  measurements  consisted  of  pure 
quartz  sand  (the  density  of  quartz  is  2.65  g/cnr  ) .  The  grain  size 
distribution  was  checked  in  standard  sieves  and  the  sand  was  found  to 
be  distributed  as  shown  in  Table  Al.  The  sand  is  very  well  sorted, 
since  approximately  70$  Qf  the  sard  grains  are  between  0.250  and 
0.177  non  in  diameter. 

Three  disk- like  samples  (diameter  5.5  in.,  thickness  «  1  in.) 
were  made  using  different  proportions  of  Scotchcast  221  and  sand.  The 
measured  and  computed  values  of  the  density  and  porosity  are  given  in 
Table  A2. 

3.  Velocity  of  Sound  in  Mixtures  of  Scotchcast  221  and  Sand 

Next  the  sound  velocity  in  each  of  the  three  samples  given 
above  was  measured.  Since  the  sound  velocity  is  fairly  difficult  to 
measure,  three  different  techniques  were  used  to  determine  the  sound 
velocity. 

METHOD  I:  The  travel  time  differences  between  a  source  and 
a  receiver  are  measured  when  the  sample  is  in  place  between  the  source 
and  receiver  and  when  it  is  removed.  The  sound  velocity  in  the  sample 
is  then  given  by 
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TABLE  A1 


DISTRIBUTION  OF  GRAIN  SIZE 
BY  WEIGHT  OF  NO.  5  SAND 


Sieve  No. 


0 

mm 


7.  by  weight 


45 

60 

80 

120 

170 

others 


0.350 
0.250  \ 
0.177 
0.125 
0.088 


5.2% 
69.87c 
18 . 87o 
4.47c 
1.87c 


TABLE  A2 

DENSITY  AND  POROSITY  VALUES 
OF  SCOTCHCAST-SAND  MIXTURE 


Mixing  Ratio  by  Weight 

Scotchcast : Sand 

Measured 

Density 

g/cm3 

Computed 

Porosity 

7. 

Pure  Scotchcast 

1.06 

1007, 

1:1 

1.51 

71.57, 

1:3 

1.86 

45.67. 
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c 


1 


) 


(A4) 


where  cwa^.er  is  the  velocity  of  sound  in  water  measured  when  the  sample 
is  removed,  d  is  the  sample  thickness,  and  At  is  the  travel  time 
difference.  While  this  method  is  not  particularly  accurate,  it  does 
have  the  ad-vantage  of  allowing  c  to  be  measured  for  a  wide  range  of 
frequencies  to  determine  dispersion.  The  velocity  of  sound  in  water 
at  72 °F  was  found  to  be 


'water 


=  1463 


The  velocity  of  sound  as  a  function  of  frequency  for  the  three 
different  samples  measured  is  given  in  Dwg.  AS- 71- 35 8 .  In  the  range 
of  frequencies  measured,  the  sound  velocity  in  pure  Scotchcast  was 
independent  of  frequency,  while  in  the  mixtures  the  velocity  appears 
to  be  slightly  frequency  dependent.  Error  analysis  of  the  measurement 
technique  indicated  that  error  for  mixture  1:1  was  about 

Ac1>^  =  ±36  m/sec  , 

while  for  mixture  1:3  the  error  was  probably 

Ac^,.  =  ±70  m/sec 

These  error  bounds  are  indicated  in  Dwg.  AS-71-358  by  the  vertical 
extensions  about  the  measured  point. 

METHOD  II:  If  the  frequency  is  varied  and  the  frequency 
difference  of  neighboring  maxima  of  the  reflected  signal  (steady  state 
part  of  a  long  pulse)  is  measured  when  the  source  and  receiver  are  on 
the  same  side  of  the  sample,  then  the  sound  velocity  can  be  calculated 
from 


105 


SOUND  VELOCITY  -  m/sec 


FREQUENCY  -  kHz 


ARL  -  UT 
AS • 71  -  358 
PJW  •  RHC 
4  -  20  •  71 


SOUND  VELOCITY  AS  A  FUNCTION  OF  FREQUENCY 
MEASURED  BY  METHOD  1 


c  =  2  dAf 


(A5) 


Of  the  three  methods  used,  this  method  should  give  the  most  accurate 
measurement  of  the  sound  velocity  if  the  absorption  is  not  too  great. 
However,  at  the  frequencies  used  here,  the  absorption  was  large  enough 
that  pulse  distortion  did  occur,  and  consequently  the  accuracy  of  the 
results  were  reduced  somewhat.  The  average  value  of  four  measurements 
for  each  of  three  samples  is  shown  in  Table  A3. 

METHOD  III:  The  time  difference  between  the  very  short 
pulses  (shorter  than  2d*c)  reflected  from  the  front  and  the  rear 
interface  of  the  sample  was  measured.  This  method  has  the  advantage 
of  measuring  the  velocity  over  twice  the  thickness  instead  of  just 
over  one  thickness,  as  Method  I.  However,  care  must  be  taken  to  allow 
for  the  phase  change  which  occurs  upon  reflection  at  the  rear  interface, 
and  the  correct  portion  of  the  two  pulse  returns  must  be  compared.  The 
average  value  of  four  measurements  for  each  of  three  samples  is  given 
in  Table  A3  under  the  heading  Method  III. 


TABLE  A3 

THE  VALUES  OF  THE  SPEED  OF  SOUND  IN 
MIXTURES  OF  SCOTCHCAST  AND  SAND 


MIXTURE 

METHOD  I 

' 

METHOD  II 

'  " 

METHOD  III 

MEAN  VALUE 

Pure  Scotch- 
cast 

1743  m/sec 

1800  m / sec 

1789  m/sec 

1777  m/sec 

1:1 

1826  m/sec 

1959  m/ sec 

1925  m / sec 

1903  m/sec 

1:3 

2253  m/sec 

2401  m/sec 

2346  m/sec 

2333  m/sec 

The  values  of  the  speed  of  sound  listed  in  Table  A3  under  Method  I  are 
the  average  values  of  the  speed  of  sound  measured  in  the  frequency 
range  100  -  500  kHz. 
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4.  Measurement  of  the  Absorption 


The  absorption  of  a  medium  is  generally  identified  as  the 
imaginary  part  of  the  complex  wave  number 

k  =  k  +  ia 
o 

Absorption  is  usually  measured  by  varying  the  source- receiver 
distance  while  maintaining  a  constant  source  output;  the  signal  loss 
after  spherical  spreading  has  been  removed  is  then  the  absorption  loss 
(assuming  a  homogeneous  medium). 

In  the  case  at  hand,  the  determination  of  the  absorption  is 
not  so  simple  since  only  one  sample  of  absorbing  material  of  fixed 
thickness  was  available.  The  absorption  can  be  determined,  however, 
by  measuring  the  reflected  level  of  a  short  pulse  and  the  front  inter¬ 
face  and  the  transmitted  level  of  this  pulse  through  the  disk.  The 
transmitted  pressure  is  given  by 

TPi  =  ,  (A6) 

where  Eq.  (A 6)  represents  the  transmitted  pressure  for  pulselengths 
Cmate ri a 1 

L  <  2d  - .  T. „  is  the  single  interface  transmission  coefficient 

c  ,  12  ° 

water 

from  water  into  the  material,  is  the  single  interface  reflection 
coefficient  from  the  material  into  water,  and  d  is  the  thickness  of 
the  sample.  Consideration  of  the  boundary  conditions  gives 


(A7) 


and 


1  +  R12  =  T12 


1  +  R21  T21 


R12  “  'R21 
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where  R  represents  the  reflection  coefficient,  and  the  subscripts 
indicate  the  same  interfaces  as  those  of  the  transmission  coefficients. 
Substituting  Eqs.  (Aj)  into  (A 6)  gives 


r  /,  „2  ^  -ad 

T-  t1  -  V  e 


The  experimental  quantities  which  have  been  measured  are  the  magnitudes 
|T|  and  |R_gi,  and  the  absolute  value  of  Eq.  (A8)  is 


|T|  -  [l  -  l^l2} 


Consequently,  a  is  given  by 


Ml1-  I*-,/]/'1-!}  . 


(A10) 


Both  |R]2|  and  |T|  have  been  measured  as  a  function  of  frequency  in  the 
range  150  to  500  kHz  for  each  of  the  three  samples.  From  these  measure¬ 
ments  a  can  be  determined  according  to  Eq.  (A10).  In  Chapter  III  of 
this  report,  it  was  pointed  out  that  if  a  has  a  first  power  frequency 
dependence,  then  k  could  be  written  as 


k  =  ko(1  +  1s)  ’ 


(All) 


where 


a'  =  xa  , 


(Al?) 


and  X  is  the  wavelength  in  the  material. 
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The  values  of  a  and  a*  for  each  of  the  mixtures  are  given  in 
Table  A.!.  The  value  of  a'  is  also  plotted  in  Dwg.  AS-71-359  as  a 
function  of  frequency.  From  the  drawing  it  appears  that  a'  is  a 
constant,  and  consequently  a  does  indeed  have  a  first- power  frequency 
dep*  ->dence. 


TABLE  A4 


THE  MEASUREMENT  OF  THE  ABSORPTION  COEFFICIENT  VERSUS  FREQUENCY 


Frequency 

i 

Pure  Scotchcast 
a  (in.-1)  ot' 

Mixture  1 
<v(in.-1) 

:  1 

a’ 

1 

Mixture  1 
Qf(in.-1) 

:  3 

a' 

15C 

! - 

0.34/in. 

I  -  " 

0. 159 

0 .39/in. 

0.194 

0 .49/in. 

0.300 

200 

0.37/in. 

0.129 

0 .38/in. 

0.142 

0  .52/in. 

0.239 

250 

0 .46/in. 

0.129 

0 .64 /in. 

0.192 

0 .45/in. 

0.165 

300 

0.65/in. 

0.152 

0 .83/in. 

0.207 

0 .68 /in. 

0.208 

350 

0.73/in. 

0.146 

0.77/in. 

0.165 

0 .84/in. 

0.220 

400 

0.81/in 

0.142 

0 .97/in. 

0.182 

0 .95 /in. 

0.218 

450 

0 .98/in. 

0 .152 

1.17/in. 

0.195 

1.07/in. 

0.218 

500 

1.29/in. 

1 

0.180 

1.46 /in. 

0.219 

1. 16 /in. 

0.213 

0 .149  ±  0.016  0  . 187  ±  0.023  0  ,223±0.034 

Mean  Mean  Mean 


If  the  mean  value  of  the  absorption  is  scaled  to  a  frequency 
of  3.5  kHz,  it  is  found  that*  the  absorption  in  dB  per  foot  for  each  of 
these  materials  is 


Pure  Scotchcast  -  0.85  dB/ft 

Mixture  1:1  1.09  dB/ft 

Mixture  1:3  1.29  dB/ft 


llu 


These  values  are  quite  typical  of  those  actually  occurring 
in  ocean  bottoms . 

5 .  Summary 

The  mean  or  best  values  of  the  data  presented  in  Tables  Ai  -  A4 
are  collected  into  Table  A5.  On  the  whole,  the  errors  in  the  measurements 
are  on  the  order  of  3$  or  less  except  for  the  absorption  coefficient,  for 
which  the  error  is  about  10$.  The  greater  dispersion  in  the  absorption 
coefficient  values  arises  because  of  the  difficulty  of  measuring  the 
reflection  coefficient  at  normal  incidence.  The  difficulty  was  simply 
geometrical;  the  precision  of  the  positioning  system  used  was  not  good 
enough  to  be  certain  that  the  transducer  beam  was  perpendicular  to  the 
sample  at  each  measurement. 

The  most  striking  feature  of  the  physical  parameters  given  in 
Table  A5  is  the  excellent  modeling  obtained.  The  values  are  quite 

typical  of  higher  speed  bottoms  ( bottom  >  Cwater^  an<i  are  8180  comPle'tel3r 
consistent  for  modeling  purposes  (that  is,  none  of  the  parameter  values 

are  anomalous  with  respect  to  the  ether  values ) . 
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